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Preface 



The purpose of this pubHcation is to describe the 
FORTRAN Hbrary subprograms and their use in either 
a FORTRAN or an assembler language program. The 
body of the publication describes the mathematical 
subprograms (which perform computations) and the 
service subprograms (which perform testing and 
utility functions). This information is intended pri- 
marily for the FORTRAN programmer; Appendix E is 
intended for the assembler language programmer. 
Additional appendixes contain algorithms ( the method 
by which a mathematical function is computed), per- 
formance statistics, descriptions of interruption and 
error procedures, storage estimates, and sample stor- 
age printouts. 

The reader should be familiar with one of the follow- 
ing publications: 

IBM System! 360 FORTRAN IV Language, Form 



IBM System/ 360 Basic FORTRAN IV Language, 
Form C28-6629 

IBM System/360 Operating System: Assembler 
Language, Form C28-6514 

IBM System/360 Model 44 Programming System: 
Assembler Language, Form C28-6811 

In addition, references are made within this publica- 
tion to information contained in the following publi- 
cations : 

IBM System/ 360 Principles of Operation, Form 
A22-6521 

IBM System/ 360 Operating System: Supervisor and 
Data Management Macro-Instructions, Form C28-6647 

IBM System/ 360 Model 44 Programming System: 
Guide to System Use, Form C28-6812 

Standard mathematical notation is used in this publi- 
cation. The reader is expected to be familiar with this 
notation and with common mathematical terminology. 



Note to users of the System/360 Operating System: The information in this pubHcation about the execution-time routines 
IHCADJST, IHCTRCH, and IHCUOPT and about the form of the program interrupt message that results from bound- 
ary alignment errors will become effective with Release 11 of the System/360 Operating System. 



Third Edition 

This publication is a major revision of Form C28-6596-1. Significant changes have been made to 
support the Model 44 Programming System FORTRAN IV Hbrary. Revisions made to text or 
tables that pertain to the use of the System/360 Operation System are indicated by vertical 
bars in the left-hand margin. 

Specifications contained herein are subject to change from time to time. Any such changes 
will be reported in subsequent revisions or Technical Newsletters. 

© International Business Machines Corporation, 1966 
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Introduction 



The FORTRAN IV library for the System/ 360 Operating 
System and the Model 44 Programming System com- 
prises two types of relocatable subprograms: mathe- 
matical subprograms and service subprograms. The 
mathematical subprograms correspond to a subpro- 
gram denned by a function statement in a Fortran 



source module. These subprograms always return one 
answer (function value) to the calling module. The 
service subprograms correspond to a subprogram de- 
fined by a subroutine statement in a Fortran source 
module. These subprograms may or may not return a 
value to the calling module. 

Calls to the library subprograms are either at the 
programmer's request or in response to program re- 
quirements. Under the System/360 Operating System, 
all calls are processed by the linkage editor, which 
takes the subprograms from the library. The library 
subprograms are then combined by the linkage editor 
with the calling module (either an object or a load 
module) into another load module which is ready for 
execution. Under the Model 44 Programming System, 
the linkage editor takes the subprograms from the 



library and combines them with the calling module into 
an executable phase. 

The library subprograms may be called in either a 
FORTRAN or an assembler language program. The next 
"ons of this publicaLion contain calling infor- 
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mation for the Fortran programmer; Appendix E con- 
tains calling information for the assembler language 
programmer. 

Subprogram Names 

The names of the subprograms described in this pub- 
lication are the same whether they are part of the 
System/ 360 Operating System or the Model 44 Pro- 
gramming System library, except for a three-character 
prefix, which is uniformly IHC for the former and 
BOA for the latter. These occupy the portion of the 
subprogram name shown as xxx in the text and •^ables. 
Thus, the subprogram here named "xxxcslog" has 
the name ihccslog under System/ 360 Operating Sys- 
tem and BOACSLOG under the Model 44 Programming 
System. 
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Mathematical Subprograms 



The mathematical subprograms suppHed in the 
FORTRAN library perform computations frequently 
needed by the applications programmer. The mathe- 
matical subprograms are called in two ways : explicitly, 
when the programmer includes the appropriate entry 
name in a source language statement (see-^ Table 1); 
and implicitly, when certain notation (e.g., raising a 
number to a power ) appears within a source language 
statement (see Table 6). 

The following text describes the individual mathe- 
matical subprograms and explains their use in a 
FORTRAN program. Detailed information about the 
actual method of computation used in each subpro- 
gram, the performance of the subprogram, interrup- 
tion and error procedures, and storage estimates can 
be found in the appendixes of this publication. 



Explicitly Called Subprograms 

Each explicitly called subprogram performs one or 
more mathematical functions. Each mathematical 
function is identified by a unique entry name that 
differs from the name of the subprogram. 

A subprogram is called whenever the appropriate 
entry name is included in a Fortran arithmetic expres- 
sion. The programmer must also supply one or more 
arguments. These arguments follow the entry name 
and are separated by commas; the list of arguments is 
enclosed in parentheses. 

For example, the source statement 

RESULT = SIN (RADIAN) 

causes the ihcsscn subprogram to be called. The sine 
of the value in radian is computed and the function 
value is stored in result. 

In the following example, the ihcssqrt subprogram 
is called to compute the square root of the value in 
amnt. The function value is then added to the value 
in STOCK and the result is stored in ans. 

ANS = STOCK + SQRT (AMNT) 

The explicitly called subprograms are described in 
the tables that make up the rest of this section. 
These tables show the general function, subprogram 
name, the Fortran library that contains the subpro- 
gram, definition, entry name(s), argument information, 
type of function value returned, and assembler re- 
quirements. The following column headings are used 
in the tables: 

General Function: This column states the nature of 
the computation performed by the subprogram. 



Subprogram Name: This column gives the module 
name of the subprogram. The name of the subprogram 
is the same whether it is in the Operating System li- 
brary or the Model 44 Programming System library 
except for the first three characters, which are uni- 
formly iHC for the forn\er and boa for the latter. These 
three characters occupy the portion of the subprogram 
name that is shown as xxx in the text and tables. 

Subset: This column indicates those subprograms 
that belong to the Fortran iv (e) library. Unless 
otherwise indicated, all subprograms that belong to 
the E library also belong to the Fortran rv ( G, H, and 
Model 44) library. 

Definition: This column gives a mathematical equa- 
tion that represents the computation. An alternate 
equation is given in those cases where there is another 
way of representing the computation in mathematical 
notation. ( For example, the square root can be repre- 
sented either as y= Vx or y = x^/^.) The definition 
for those subprograms that accept complex arguments 
contains the notation z = Xi 4- X2i. 

Entry Name: This column gives the entry name that 
the programmer must use to call the subprogram. A 
subprogram may have more than one entry name; the 
particular entry name used depends upon the compu- 
tation to be performed. For example, the ihcsscn sub- 
program has two entry names: cos and sin. If the 
cosine is to be computed, entry name cos is used; if 
the sine is to be computed, entry name sin is used. 
Argument Number: This column gives the number 
of arguments that the programmer must supply. 

Argument Type: This column describes the mode 
and length of the argument, integer, real, and 
COMPLEX represent the type of number; the notation 
*4, *8, and *16 represent the size of the argument in 
storage locations. 

Note: In Fortran rv (e), a real argument corre- 
sponds to the real *4 argument, and a double- 
precision argument corresponds to the real *8 argu- 
ment. Complex arguments cannot be used with a 
FORTRAN IV (e) compiler. 

Argument Range: This column gives the valid range 
for an argument. If the argument is not within this 
range, an error message is issued and execution of this 
load module is terminated. Appendix C contains a 
description of the error messages. 

Function Value Type: This column describes the 
type of function value returned by the subprogram. 
The notation used is the same as that used for the 
argument type. 



Table 1. Explicitly Called Mathematical Subprograms 



General Function 


Specific Function 


Subprogram Name 


Entry Name(s) 


Logarithmic and exponential 
subprograms ( described in Table 2 ) 


Common and natural logarithm 


xxxCLLOG* 
xxxCSLOG* 
xxxLLOG 
xxxSLOG 


CDLOG 
CLOG 

DLOG, DLOGIO 
ALOG, ALOGIO 


Exponential 


xxxCLEXP* 
xxxCSEXP* 
xxxLEXP 
xxxSEXP 


CDEXP 
CEXP 
DEXP 
EXP 


Square root 


xxxCLSQT 

XXX\jOO\/ 1 

xxxLSQRT 
xxxSSQRT 


CDSQRT 

DSQRT 
SQRT 


Trigonometric subprograms 
( described in Table 3 ) 


Arcsine and arccosine 


xxxLASCN* 
xxxSASCN* 


DARSIN, DARCOS 
ARSIN, ARCOS 


Arctangent 


IHCLATAN 
XXXLATN2* 
IHCSATAN 
XXXSATN2* 


DATAN 

DATAN, DATAN2 
ATAN 
ATAN, ATAN2 


Sine and cosine 


xxxCLSCN* 
xxxCSSCN* 
xxxLSCN 
xxxSSCN 


CDSIN, CDCOS 
CSIN, CCOS 
DSIN, DCOS 
SIN, COS 


Tangent and cotangent 


xxxLTNCT* 
xxxSTNCT* 


DTAN, DCOTAN 
TAN, COTAN 


Hyperbolic function subprograms 
( described in Table 4 ) 


Hyperbolic sine and cosine 


xxxLSCNH* 
xxxSSCNH* 


DSINH, DCOSH 
SINH, COSH 


Hyperbohc tangent 


xxxLTANH 
xxxSTANH 


DTANH 
TANK 


Miscellaneous subprogiams 
( described in Table 5 ) 


Absolute value 


xxxCLABS* 
xxxCSABS* 


CDABS 
CABS 


Error function 


xxxLERF* 
xxxSERF* 


DERF, DERFC 
ERF, ERFC 


Gamma and log-gamma 


xxxLGAMA* 
xxxSGAMA* 


DCAMMA, DLGAMA 
GAMMA, ALGAMA 


Maximum and minimum value 


xxxFMAXD 

xxxFMAXI 

xxxFMAXR 


DMAXl, DMINl 

AMAXO, AMINO, MAXO, MINO 

AMAXl, AMINl, MAXl, MINI 


Modular arithmetic 


IHCFMODI 
IHCFMODR 


MOD 

AMOD, DMOD 


Truncation 


IHCFAINT 
IHCFIFIX 


AINT 
IDINT, INT 


*Not available in FORTRAN IV (E) 



Assembler Requirements: This column gives the reg- 
isters used by the subprogram and the minimum save 
area that the assembler language programmer must 
supply. For example, the assembler requirements for 
the a:xxcssQT subprogram are: 
registers 0, 2(4) 
save area 9F 
This information specifies that: 

1. The function value is found in floating-point reg- 
isters and 2. 



2. Floating-point register 4 is used for intermediate 
computation. 

3. The save area must be at least nine full-words in 
length. 

Detailed information for the assembler language pro- 
grammer is given in Appendix E. 

Note: In the following tables, the approximate value 
of 2'^ • TT is .82354966406249996D+06; the approximate 
value of 2P'> • tt is .35371188737802239D + 16. 
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• Table 2. Lxjgarithmic and Exponential Subprograms 





General 


.SubrirnPTam 


Sub- 
set 


Definition 


Entry 
Name 


Argument(s) 


Function 
Value 
Type' 


Assembler 
Require- 
ments 




Function Name 


No. 


Type' 


Range 




Common 
and natural 
logarithm 


xxxCLLOG 


No 


y = pvlogc (z) 
See Note 2 


CDLOG 


1 


complex 
*16 


z ^ -1- Oi 

See Note 3 


complex 
*16 


registers 
0,2 

save area 
8F 




xxxCSLOG 


No 


y=pvlogc (z) 
See Note 2 


CLOG 


1 


complex 
*8 


z ^ -f Oi 

See Note 3 


complex 
*8 


registers 
0,2 

save area 
8F 




xxxLLOG 


Yes 


y=logeXor 
y=lnx 


DLOG 


1 


real *8 


x>0 


real *8 


registers 
0(2) 
save area 
9F 




y=logiox 


DLOGIO 


1 


real *8 


x>0 


real *8 


registers 
0(2) 
save area 
9F 




xxxSLOG 


Yes 


y=logeX or 
y=lnx 


ALOG 


1 


real *4 


x>0 


real *4 


registers 
0(2) 
save area 
5F 




y=logiox 


ALOGIO 


1 


real *4 


x>0 


real *4 


registers 
0(2) 
save area 
5F 


1 


Exponential 


xxxCLEXP 


No 


y=e^ 

See Note 4 


CDEXP 


1 


complex 
*16 


xi < 174.673 
N<(2-.,r) 


complex 
*16 


registers 
0,2 

save area 
8F 


1 


xxxCSEXP 


No 


y=e^ 

See Note 4 


CEXP 


1 


complex 
*8 


xi < 174.673 

|x2|<(2=»-7r) 


complex 
*8 


registers 
0,2 

save area 
8F 




xxxLEXP 


Yes 


y=e^ 


DEXP 


1 


real *8 


X ^ 174.673 


real *8 


registers 
0(2) 
save area 
9F 




xxxSEXP 


Yes 


y=e^ 


EXP 


1 


real *4 


X ^ 174.673 


real *4 


register 
save area 
12F 




Square root 


xxxCLSQT 


No 


y= Vz or 

y=z^/^ 


CDSQRT 


1 


complex 
*16 


any complex 
argument 
See Note 3 


complex 
*16 


registers 
0,2(4) 
save area 
9F 




xxxCSSQT 


No 


y= Vz or 


CSQRT 


1 


complex 
*8 


any complex 
argument 
See Note 3 


complex 
*8 


registers 
0,2(4) 
save area 
9F 




xxxLSQRT 


Yes 


y= Vxor 


DSQRT 


1 


real ^8 


x^O 


real *8 


registers 
0(2,4) 
save area 
5F 




xxxSSQRT 


Yes 


y= Vx or 


SQRT 


1 


real *4 


x^O 


real *4 


registers 
0(4) 
save area 
5F 


1 


Notes: 

1. In FORTRAN IV (E), a real argument corresponds to the REAL *4 argument, and a double-precision argument corresponds to 
the REAL *8 argument. 

2. PV = principal value. The answer given is from that point where the imaginary part (yi-) lies between — tt and +Tr. More specifi- 
cally: — TT < ya < IT, unless xi < and X2 = —0, in which case, y2 = — tt. 

3. Floating-point overflow can occur. 

4. Where z is a complex number of the form Xi-t-xai. 



• Table 3. Trigonometric Subprograms 





General 
Function 


ijubprogram 
Name 


Sub- 
set 


Definition 


Entry 
Xame 


Argument(H) 


Function 
Value 
Type' 


Assembler 
Require- 
ments 




No. 


Type' 


Range 




Arcsine and 
arccosine 


xxxLASCX 


No 


y = arcsin (x) 


DARSIN 


1 


real *8 


1x1 ^1 


real 
*8(in 
radians ) 


registers 
(2,4) 
save area 
13F 




y = arccos (x) 


DARCOS 


1 


real *8 


Ixl^l 


real 
*8(in 
radians ) 


registers 
0(2,4) 
save area 
13F 




xxxSASCX 


Xo 


y = arcsin (x) 


ARSIN 


1 


real *4 


Ixl^l 


real 

*4(in 
radians ) 


registers 
0(2,4) 
save area 
lOF 




y = arccos (x) 


ARCOS 1 


real *4 


ixi^l 


real 
*4 (in 
radians) 


registers 
0(2,4) 
save area 
lOF 




Arctangent 


IHCLATAN 


See 

Note 

2 


y = arctan (x) 


DATAN 1 1 real *8 i any real 
i i argmnent 

i 1 ^ 


real 
*8(in 
radians ) 


registers 
(2,4,6) 
save area 
5F 




xxxLATX2 


See 

Xote 

2 


y = arctan (x) 


DATAN 


1 i real *8 

i 
i 


any real 

argimient 


real 
*8(in 
radians ) 


registers 
0(2,4,6) 
save area 
5F 




i 


/xA j DATAX2i 2 
y = arctan (-) | | 

t 1 


real *8 


any real real 
arguments *8 (in 
( except 0, ) , radians ) 


registers 
0(2,4,6) 
save area 
5F 




IHCSATAN 


See 

Xote 

2 


y = arctan (x) 


ATAN 


1 


real *4 


any real 
argument 


real 
*4(in 
radians ) 


registers 
0(2,4,6) 
save area 
5F 




XXXSATN2 


See 

Xote 
2 


y = arctan ( x ) 


ATAN 


1 


real *4 


any real 
argument 


real 
*4 (in 
radians ) 


registers 
0(2,4,6) 
save area 
5F 




y = arctan I — 1 


ATAN2 


2 


real *4 


any real 
arguments 

(except 0,0) 


real 
*4 (in 

radians ) 


registers 
0(2,4,6) 
save area 
5F 


1 


Sine and 
cosine 


xxxCLSCN 


Xo 


y = sin (z) 
See Note 4 


CDSIN 


1 


complex 
*16 (in 
radians ) 


!x,|<(2'^-7r) 

|x.l^ 174.673 


complex 
*16 


registers 
0,2(4) 
save area 
9F 


1 


y = cos ( z ) 
See Xote 4 


CDCOS 


1 


complex 
*16 (in 
radians ) 


|xi|<(2«'-7r) 
jx^i ^ 174.673 


complex 
*16 


registers 
0,2(4) 
save area 
9F 


1 


xxxCSSCN 


No 


y = sin (z) 
See Xote 4 


CSIN 


1 


complex 
*8(in 
radians ) 


!xJ ^ 174.673 


complex 
*8 


registers 
0,2(4) 
save area 
9F 


1 


y = cos (z) 
See Xote 4 


CCOS 


1 


complex 
*8(in 

radians ) 


'x>!<(2"-7r) 

\xJ ^ 174.673 


complex 
*8 


registers 
0,2(4) 
save area 
9F 
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Table 3. Trigonometric Subpr 


jgrams 


( Continued ) 














General 
Function 


Subprogram 
A'ame 


Sub- 
set 


Definition 


Entry 
Name 


Argument(s) 


Function 
Value 
Type' 


Assembler 
Require- 
ments 


No. 


Type^ 


Range 


Sine and 

cosine 

( continued ) 


xxxLSCN 


Yes 


y = sin (x) 


DSIN 


1 


real 
*8 (in 
radians ) 


|x|<(2^-7r) 


real *8 


registers 
0(2,4) 
save area 
5F 


y = cos (x) 


DCOS 


1 


real 
*8(in 
radians ) 


|x|<(2«'.7r) 


real *8 


registers 
0(2,4) 
save area 
5F 


xxxSSCN 


Yes 


y = sin (x) 


SIN 


1 


real 

*4 (in 
radians ) 


!xi<(2--,r) 


real *4 


registers 
0(2,4) 
save area 
5F 


y = cos (x) 


COS 


1 


real 
=^4 (in 
radians ) 


ix|<(2'^-:r) 


real *4 


registers 
0(2,4) 
save area 
5F 


Tangent and 
cotangent 


xxxLTNCT 


No 


y = tan (x) 


DTAN 


1 


real 
*8(in 
radians ) 


ix|<(2-.-7r) 
See Note 3 


real *8 


registers 
0(2,4,6) 
save area 
5F 


y = cotan (x) 


DCOTAN 


1 


real 
*8(in 
radians ) 


Ixl<(2«'-7r) 

See Note 3 


real *8 


registers 
0(2,4,6) 
save area 
5F 


xxxSTNCT 


No 


y = tan (x) 


TAN 


1 


real 
*4(in 
radians ) 


|x|<(2^«-7r) 
See Note 3 


real *4 


registers 
0(2,4) 
save area 
5F 


y = cotan ( x ) 


COTAN 


1 


real 
*4(in 
radians ) 


|x|<(2^^-7r) 

See Note 3 


real *4 


registers 
0(2,4) 
save area 
5F 


Notes: 

1. In FORTRAN IV (E), a real argument corresponds to the REAL *4 argument, and a double-precision argument corresponds to 
the REAL *8 argument. 

2. Instead of the IHCLATAN and IHCSATAN subprograms contained in the FORTRAN IV (E) library, the FORTRAN IV library 
contains the xxxLATN2 and xxxSATN2 subprograms. 

3. The argument for the cotangent functions may not be near a multiple of rr; the argument for the tangent functions may not be near 
an odd multiple of 7r/2. 

4. Where z is a complex number of the form Xi + xai. 
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Table 4. Hyperbolic Function Subprograms 



General 
Function 


Subprogram 
Name 


Sub- 
set 


Definition 


Entry 
Name 


Argument(s) 


Function Assembler 
Value Require- 
Type- ments 


No. 


Type- 


Range 


Hyperbolic 
sine and 
cosine 


xxxLSCNH 


No 


e'' — e"'' 

y- 2 


DSINH 


1 


real *8 


|xj<174.673 


real *8 registers 
0(2,4) 
save area 
9F 


e"" -1- e"'" 
y- 2 


DCOSH 


1 


real *8 


[xj<174.673 


real *8 registers 
0(2,4) 
save area 
9F 


xxxSSCNH 


No 


e" — e"'" 


SINH 


1 


real *4 


|xi<174.673 


real *4 registers 
0(2,4) 
save area 
9F 


^' 2 


e^-l-e"'' 
y- 2 


COSH 


1 


real *4 


|x|<l74.673 


real *4 registers 
0(2,4) 
save area 
9F 


Hyperbolic 
tangent 


xxxLTANH 


Yes 


e"" — e"" 
^ e'^ + e-'' 


DTANH 


1 


real *8 


any real 
argument 


real *8 registers 
0(2) 
save area 
5F 


xxxSTANH 


Yes 


e" — e"" 


TANH 


1 


real *4 


any real 
argument 


real *4 registers 
0(2) 
save area 
5F 


^ e'^ + e"" 


Note: 

1. In FORTRAN IV (e), a real argument corresponds to the REAL *4 argument, and a double-precision argument corresponds 
to the REAL *8 argument. 



rauie 5. Miscellaneous Mathematical Subprograms 



General 
Function 


Subprogram 
Name 


Sub- 
set 


Definition 


Entry 
Name 


Argument(s) 


Function Assembler 
Value Require- 
Type^ ments 


No. 


Type^ 


Range 


- 
Absolute 
value 

■ 


xxxCLABS 


No 


y=izj = (xi' + X2')''' 


CDABS 


1 


complex 
*16 


any complex 
argument ' 
See Note 2 


real *8 registers 
0,2(4) 
save area 
8F 


xxxCSABS 


No 


y=|z| = (xiHx/)^^= 


CABS 


1 


complex 
*8 


any complex 
argument 
See Note 2 


real *4 registers 
0,2(4) 
save area 
8F 


Error 
function 


vxvLERF 


No 


y— — 1 e " du 


DERF 


1 


real *8 


any real 
argument 


real *8 registers 
0(2,4,6) 
save area 
IIF 


2 fV. 
v= — = 1 e du 

y=l — erf (x) 


DERFC 


1 


real *8 


any real 
argument 


real *8 registers 
0(2,4,6) 
save area 
IIF 
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Table 5. Miscellaneous Mathematical Subprt)granis (Continued; 



General 
Function 


Siibprojfram 
Name 


Sub- 
set 




Entry 


Arguinent(s) 


Function 
Value 
Type^ 


Assembler 
Require- 
ments 


Definition I Name 


No. \ Type' 


Range 


Error 
function 
( continued ) 


xxxSERF 


No 


2 /•" -" ; ERF ; 1 ; real *4 

V— — 1 e " du ; i 

i 1 i 


any real 
argument 


real *4 


registers 
0(2,4,6) 
save area 
IIF 


2 /•== . JeRFC 
y— — 1 e " du 

VttJ ^ 

y=l — erf (x) j 


1 


real *4 


any real 
argument 


real *4 


registers 
0(2,4,6) 
save area 
IIF 


Gamma and 
log-gamma 


xxxLGAMA 


No 


r « DGAMMA 
y= 1 u'^'^e"" du 

1 


1 


real *8 


X > 2-=^-^^' and 
X < 57.5744 


real *8 


registers 
0(2,4,6) 
save area 
IIF 


y = loge r (x) or 
y=loge 1 u'^'^e"" du 


DLGAMA 


1 


real *8 


X > and 

x< 4.2913 -10" 


real *8 


registers 
0(2,4,6) 
save area 
IIF 


xxxSGAMA 


No 


y = 1 u" ' e " du 


GAMMA 


1 


real *4 


X > 2-'^^= and 
X < 57.5744 


real *4 


registers 
0(2,4,6) 
save area 
IIF 


y = log<. r (x) or 
y=loge 1 u'' ^ e " du 


ALGAMA 


1 


real *4 


X > and 

x< 4.2913- 10'' 


real *4 


registers 
0(2,4,6) 
save area 
IIF 


Maximum 
and 

minimum 
values 


xxxFMAXD 


Yes 


y = max (xi , . . . , Xn) 


DMAXl 


>2 


real *8 


any real 
arguments 


real *8 


register 
save area 
9F 


y = min (xi, . . . ,Xn) 


DMINl 


^2 


real *8 


any real 
arguments 


real *8 


register 
save area 
9F 


xxxFMAXI 


Yes 


y=max (xi , . . . ,xn) 


AMAXO 


^2 


integer 

*2^ *4 


any integer 
arguments 


real *4 


register 
save area 
9F 


MAXO 


^2 


integer 
*2,*4 


any integer 
arguments 


integer 
*2,*4 


register 
See Note 3 
save area 
9F 


y = min (xi , . . . , Xn) 


AMINO 


^2 


integer 

*2,*4 


any integer 
arguments 


real *4 


register 
save area 
9F 


MINO 


^2 


integer 
*2^ *4 




any integer 
arguments 


integer 
*2, *4 


register 
See Note 3 
save area 
9F 
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• Table 5. Miscellaneous Mathematical Subprograms ( Continued ) 



General 
Function 


Subprogram 
Name 


Sub- 
set 


Definition 


Entry 
Name 


Argument(s) 


Function 
Value 
Type' 


'- 
Assembler 
Require- 
ments 


No. Type' 


Range 


Maximum 

and 

Minimum 

Values 

( continued ) 


xxxFMAXR 


Yes 


y=max (xi , . . . , Xn) 


AMAXl 


^2 real*4 


any real 
arguments 


real *4 


register 
save area 
9F 


MAXl 


^2 real *4 


any real 
arguments 


integer 
*2, *4 


register 
See Note 3 
save area 
9F 


y=min (xi , , . T,Xn) 


AMINl 


^2 real*4 


any real 
arguments 


real *4 


register 
save area 
9F 


MINI 


^2 real*4 


any real 
arguments 


integer 
*2,*4 


register 
See Note 3 
save area 
9F 


Modular 
arithmetic 


IHCFMODI 


See 

Note 

4 


y = Xi ( modulo X2 ) 
See Note 5 


MOD 


2 integer 


X2 7^0 

See Note 6 


integer 
*2,*4 


register 
See Note 3 
save area 
9F 


IHCFMODR 


See 

Note 

4 


y = xi (modulo X2) 
See Note 5 


AMOD 


2 real *4 


X2^0 

See Note 6 


real *4 


register 
save area 
9F 


y = xi (modulo X2) 
See Note 5 


DMOD 


2 real *8 


X2 7^0 

See Note 6 


real *8 


register 
save area 
9F 


Truncation 


IHCFAINT 


See 

Note 

4 


y = (signx) 'r 
where n is the largest 
integer ^ |x| 


AINT 


I real *4 


any real 
argument 


real *4 


register 
save area 
9F 


IHCFIFIX 


See 

Note 

4 


y = (signx) -n 
where n is the largest 
integer ^ |xj 


IDINT 


1 real *8 


any real 
argument 


integer 


register 
See Note 3 
save area 
9F 


INT 


1 real *4 


any real 
argument 


integer 


register 
See Note 3 
save area 
9F 


Notes: 

1. In FORTRAN IV (E), a real argument corresponds to the REAL *4 argument, and a double-precision argument corresponds 
to the REAL *8 argument. 

2. Floating-point overflow can occur. 

3. The result is stored in general register 0. 

4. The coding that performs this function is out-of-line in FORTRAN IV (E) and in-line in FORTRAN IV. Out-of-line coding is 
taken from the FORTRAN library by the linkage editor and processed with the calhng module. In-line coding is inserted by the 
FORTRAN compiler at the point in the source module where the function is referenced. This means that the in-line functions 
are available in FORTRAN IV by using the appropriate entry name but that they are not part of the library. 

5. The expression xi (modulo xs) is defined as x, - — 1 • X2, where the brackets indicate that an integer is used. The largest in- 

L X2 J 

teger whose magnitude does not exceed the magnitude of — is used. The sign of the integer is the same as the sign of — . 

X2 X2 

6. If X2 = 0, then the modulus function is mathematically undefined. In addition, a divide exception is recognized and an interrup- 
tion occurs. (A detailed description of the interruption procedure is given in Appendix C.) 
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Implicitly Called Subprograms 

The implicitly called subprograms perform operations 
required by the appearance of certain notation in a 
FORTRAN source statement. When a number is to be 
raised to a power or when multiplication and division 
of complex numbers are to be performed, the Fortran 
compiler generates the instructions necessary to call 
the appropriate subprogram. For example, if the fol- 
lowing source statement appears in a source module, 

ANS = BASE**EXPON 

where base and expon are values of the form real *4, 
the xxjcfrxpr subprogram is called by the Fortran 
compiler. 

The implicitly called subprograms in the Fortran 
library are described in Table 6. This table shows the 



general function, subprogram name, the Fortran 
library that contains the subprogram, implicit function 
reference, entry name, argument information, type of 
function value returned, and assembler requirements. 
The column headed "Implicit Function Reference" 
gives a sample source statement that might appear in 
a FORTRAN source module and cause the subprogram 
to be called. The rest of the column headings in Table 
6 have the same meaning as those used with the ex- 
plicitly called subprograms. 

The action taken within the subprogram depends 
upon the type of base and exponent used. Tables 7 
through 10 show the result of an exponentiation per- 
formed with the different combinations and values of 
base and exponent. In these tables, I and J are integers; 
A and B are real numbers; C is a complex number. 



I Table 6. Implicitly Called Mathematical Subprograms 



General 
Function 


Subprogram 
Name 


Sub- 
set 


Implicit 
Fvmction 
Reference^ 


Entry= 
Name 


Argument(s) 


Function 
Value 
Type« 


Assembler 
Requirements 


No. 


Type« 


Multiply and 
divide complex 
numbers 


xxxCLAS 


No 


y = zi* z-2 


CDMPY# 


2 


complex *16 


complex *16 


registers 0, 2 (4, 6) 
save area 5F 


y = Zi/z2 


CDDVD# 


2 


complex *16 


complex *16 


registers 0, 2 (4, 6) 
save area 5F 


xxxCSAS 


No 


y = Zi* Z2 


CMPY# 


2 


complex *8 


complex *8 


registers 0, 2 (4, 6) 
save area 5F 


y = Zi/z2 


CDVD# 


2 


complex *8 


complex *8 


registers 0, 2 (4, 6) 
save area 5F 


Raise an integer 
to an integral 
power 


xxxFIXPI 


Yes 


y = i**] 


FIXPI# 


2 


i = integer *4 
j = integer *4 


integer *4 


register 
See Note 4 
save area 18F 


Raise a real 
number to an 
integral power 


xxxFRXPI 


Yes 


y = a**] 


FRXPI# 


2 


a = real *4 
j = integer 


real *4 


register 
save area 18F 


xxxFDXPI 


Yes 


y = a**j 


FDXPI# 


2 


a = real *8 
j = integer* 4 


real *8 


register 
save area 18F 


Raise a real 
number to a 
real power 


xxxFPXPR 


Yes 


y = a**b 


FRXPR# 


2 


a = real *4 
b = real *4 


real *4 


register 
save area 18F 


xxxFDXPD 


Yes 


y = a**b 


FDXPD# 


2 


a = real *8 
b = real ^8 


real *8 


register 
save area 18F 


Raise a complex 
number to an 
integral power 


xxxFCDXI 


No 


y = z**] 


FCDXI# 


2 


z = complex 

*16 
j = integer 


complex *16 


register 
save area 18F 


xxxFCXPI 


No 


y = z**j 


FCXPI# 


2 


z = complex 

*8 
j = integer 


complex *8 


register 
save area 18F 


Notes: 

1. This is only a representation of a FORTRAN statement; it is not the only way the subprogram may be called. 

2. This name must be used in an assembler language program to call the subprogram; the character # is a part of the name and 
must be included. 

3. In FORTRAN IV (E), a real argument corresponds to the REAL *4 argument and a double precision argument corresponds to 
the REAL *8 argument. 

4. The result is stored in general register 0. 
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Table 7. Exponentiation With Integer Base and Exponent 



Table 9. Exponentiation with Real Base and Exponent 



Uusc \1.J 


Exponent ( J ) j 


J>0 


J = o 


J<0 


I>0 


Compute the 


Function 


Function value = 1 




function value 


value = 1 


if 1 = 1 
Otherwise, function 
value = 


1 = 


Function 


Error message 


Error message 




value = 


IHC241I; 


IHC241I; 






or OA24iI 


or OA241I 




Compute the 


Function 


Function value = — 1 


KO 


function value 


value = 1 


ifl= -land 

if J is an odd 

number 
Function value = 1 

ifl= -land 

if J is an even 

number 
Otherwise, function 

value = 



Table 8. Exponentiation with Real Base and Integer Exponent 



Base (A) 


Exponent ( J ) 


J>0 


J = o 


J<0 


A>0 


Compute the 
function value 


Function 
value = 1 


Compute the 
function value 


A = 


Function 
value = 


Error message 
IHC242I or 
IHC2431; 
or OA242I 
or OA243I 


Error message 
IHC242I or 
IHC243I; 
or OA242I 
or OA243I 


A<0 


Compute the 
function value 


Function 
value = 1 


Compute the 
function value 



Base (A) 


Exponent (B) 


B>0 


B = 


B<0 


A>0 


Compute the 
function value 


Function 
value = 1 


Compute the 
function value 


A = 


Function 
value = 


Error message 
IHC244I or 
IHC245I; 
or OA244I 
or OA245I 


Error message 
IHC244I or 
IHC245I; 
or OA244I 
or OA245I 


A<0 


Error message 
IHC253I or 
IKC263I 


Function 
value = 1 


Error message 
IHC253I or 
IHC263I 



Table 10. Exponentiation with Complex Base and Integer 
Exponent 



Base (C) 
C = R + Ri 




Exponent (J) 




J>0 


J = o 


J<0 


R > and 
Ri>0 


Compute the 
function value 


Function 
value =1 + Oi 


Compute the 
function value 


R > and 
Ri = 


Compute the 
function value 


Function 
value =1 + Oi 


Compute the 
function value 


R > and 
Ri<0 


Compute the 
function value 


Function 
value =1 + Oi 


Compute the 
function value 


R = and 
Ri>0 


Compute the 
function value 


Function 
value = 1 + Oi 


Compute the 
function value 


R = and 
Ri = 


Function value 
+ Oi 


Error message 
IHC246I or 
IHC247I; 
or OA246I 
or OA247I 


Error message 
IHC246I or 
IHC247I; 
or OA246I 
or OA247I 


R = and 
Ri<0 


Compute the 
function value 


Function 
value =1 + Oi 


Compute the 
function value 


R < and 
Ri>0 


Compute the 
function value 


Function 
value = 1 + 01 


Compute the 
function value 


R < and 
Ri = 


Compute the 
function value 


Function 
value =1 + Oi 


Compute the 
function value 


R < and 
Ri<0 


Compute the 
function value 


Function 
value = 1 + Oi 


Compute the 
function value 
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Service Subprograms 



The service subprograms supplied in the Fortran 
library are divided into two groups: one group tests 
machine indicators and the other group performs 
utility functions. Service subprograms are called by 
using the appropriate entry name in a Fortran source 
language call statement. 



Machine Indicator Test Subprograms 

The machine indicator subprograms (xxxfslit, 
xxxFOVER, and xxxfdvch) test the status of pseudo in- 
dicators and may retinrn a value to the calling program. 
When the indicator is zero, it is off; when the indica- 
tor is other than zero, it is on. In the following descrip- 
tions of the subprograms, i represents an integer 
expression and / represents an integer variable. 

xxxFSLIT Subprogram 

The XXXFSLIT subprogram is used to alter, test, and/ or 
record the status of pseudo sense lights. Either of two 
entry names (slite or slitet) is used to call the sub- 
program. The particular entry name used in the call 
statement depends upon the operation to be per- 
formed. 

If the four sense lights are to be turned off or one 
sense light is to be turned on, entry name slite is used. 
The source language statement is: 
CALL SLITE(i) 
where i has a value of 0, 1, 2, 3, or 4. 

If the value of i is 0, the four sense lights are turned 
off; if the value of i is 1, 2, 3, or 4, the corresponding 
sense light is turned on. If the value of i is not 0, 1, 2, 
3, or 4, an error message is issued and execution of 
this module (or phase) is terminated. (This error 
message is explained in Appendix C. ) 

If a sense light is to be tested and its status recorded, 
entry name slitet is used. The source language state- 
ment is : 

CALL SLITET (i,/) 
where: 

/ has a value of 1, 2, 3, or 4, and indicates which 
sense light to test. 

/ is set to 1 if the sense light was on; or to 2 if the 
sense light was off. 

If the value of i is not 1, 2, 3, or 4, an error message 
is issued and execution of this module (or phase) is 
terminated. ( This error message is explained in Appen- 
dix C.) 



xxxFOVER Subprogram 

The XXXFOVER subprogram tests for an exponent over- 
flow or underflow exception and returns a value that 
indicates the existing condition. After testing, the over- 
flow indicator is turned off. This subprogram is called 
by using the entry name overfl in a call statement. 
The source language statement is: 

CALL OVERFL (/) 
where: 

/ is set to 1 if a floating-point overflow condition 
exists; to 2 if no overflow or underflow condition 
exists; or to 3 if a floating-point underflow condi- 
tion exists. A detailed description of each exception 
is given in Appendix C. 

xxxFDVCH Subprogram 

The XXXFDVCH subprogram tests for a divide-check 
exception and returns a value that indicates the exist- 
ing condition. After testing, the divide-check indicator 
is turned off. This subprogram is called by using entry 
name dvchk in a call statement. The source language 
statement is: 

CALL DVCHK (/) 
where: 

/ is set to 1 if the divide-check indicator was on; or 
to 2 if the indicator was off. A detailed description of 
the divide-check exception is given in Appendix C. 

Utility Subprograms 

The utility subprograms perform two operations for 
the FORTRAN programmer: they either terminate execu- 
tion (xxxFExn) or dump a specified area of storage 

( XXXFDUMP ) . 

xxxFEXIT Subprogram 

tion ( XXXFEXIT subprogram terminates execution of this 
load module (or phase) and returns control to the 
operating system. (This subprogram performs a func- 
tion similar to that performed by the stop statement. ) 
The xxxFExrr subprogram is called by using the entry 
name exit in a call statement. The source language 
statement is: 

CALL EXIT 

xxxFDUMP Subprogram 

The XXXFDUMP subprogram dumps a specified area of 
storage. Either of two entry names (dump or pdump) 
can be used to call the subprogram. The entry name 
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is followed by the limits of the area to be dumped and 
the format specification. The entry name used in the 
CALL statement depends upon the nature of the dump 
to be taken. 

If execution of this load module ( or phase ) is to be 
terminated after the dump is taken, entry name dump 
is used. The source language statement is: 

CALL DUMP (fli, h, fi,..., an, K, U) 
where: 

a and h are variables that indicate the limits of stor- 
age to be dumped ( either a or h may represent the 
upper or lower limits of storage). 
/ indicates the dump format and may be one of the 
integers given in Table 11. The formats available 
depend upon the compiler in use. A sample printout 
for each format is given in Appendix F. 

Table 11. The xxxFDUMP Subprogram Format Specifications 



FORTRAN IV (E) 


FORTRAN IV 


specifies hexadecimal 

4 specifies INTEGER 

5 specifies REAL 

6 specifies DOUBLE 
PRECISION 


specifies hexademical 

1 specifies LOGICAL *1 

2 specifies LOGICAL *A 

3 specifies INTEGER *2 

4 specifies INTEGER *4 

5 specifies REAL *4 

6 specifies REAL *8 

7 specifies COMPLEX *8 

8 specifies COMPLEX *16 

9 specifies hteral 



If execution is to be resumed after the dump is taken, 
entry name pdump is used. The source language state- 
ment is : 

CALL PDUMP (ci, K /i, . . . , fln, hn, U) 

where a, b, and / have the same meaning as explained 
previously. 

Programming Considerations 

A load module (or, in the Model 44 Programming 
System, a member of the phase library) may occupy 
a diflFerent area of storage each time it is executed. To 
ensure that the appropriate areas of storage are dumped, 
the following conventions should be observed. 

Note: In the following examples, a is a variable in 
COMMON, B is a real number, and the array table is 
dimensioned as: 

DIMENSION TABLE (20) 

If an array and a variable are to be dumped at the 
same time, a separate set of arguments should be 
used for the array and for the variable. The specifica- 



tion of limits for the array should be from the first 
element in the array to the last element. For example, 
the following call to the ihcfdump subprogram could 
be used to dump table and b in hexadecimal format 
and terminate execution after the dump is taken: 

CALL DUMP (TABLE (1), TABLE (20), 0, B, B, 0) 

If an area of storage in common is to be dumped 
at the same time as an area of storage not in common, 
the arguments for the area in common should be given 
separately. For example, the following call to the 

-..^-vr^^^r^ ,x, „„l, ^^ lA U^ ,,-,^,1 4-^ A,.^^ 4-U^ 

AAArUKJjM.r aU.I.JI^l.\Jii.i.Ciill V^VJUIU. UC LUC»a HJ U.L1111L; LilO 

variables a and b in real *8 format without terminat- 
ing execution: 

CALL PDUMP ( A,A,6,B,B,6 ) 

If variables not in common are to be dumped, each 
variable must be listed separately in the argument 
list. For example, if r, p, and q are defined implicitly 
in the program, the statement 

CALL PDUMP (R,R,5,P,P,5,Q,Q,5) 

should be used to dump the three variables. If the 
statement 

CALL PDUMP (R,Q,5) 

is used, all main storage between r and Q is dumped, 
which may or may not include p, and may include 
other variables. 

If an array and a variable are passed to a subroutine 
as arguments, the arguments in the call to the xxxfdump 
subprogram in the subroutine should specify the 
parameters used in the definition of the subroutine. 
For example, if the subroutine subi is defined as: 

SUBROUTINE SUBI(X,Y) 
DIMENSION X(10) 

and the call to subi within the source module is : 

DIMENSION A (10) 



CALL SUBI (A, B) 

then the following statement in the subroutine should 
be used to dump the variables in hexadecimal format 
without terminating execution: 

CALL PDUMP (X(1),X(10),0,Y,Y,0) 

If the statement 

CALL PDUMP (X(1),Y, 0) 

is used, all storage between a(i) and y is dumped, 
due to the method of transmitting arguments. 
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Appendix A. Algorithms 



Appendix A contains information about the computations used in the explicitly 
called mathematical subprograms. This information is arranged in alphabetical 
order, according to the module (or phase) name of the subprogram. The entry 
names associated with each subprogram are given in parentheses after the module 
(or phase) name. 

The information for each subprogram is divided into two parts. The first part 
describes the algorithm used; the second part describes the effect of an argument 
error upon the accuracy of the answer returned. 

The presentation of each algorithm is divided into its major computational steps; 
the formulas necessary for each step are supplied. Some of the formulas are 
widely known; those that are not so widely known are derived from more common 
formulas. The process leading from the common formula to the computational 
formula is sketched in enough detail so that the derivation can be reconstructed 
by any one who has an understanding of higher mathematics and access to the 
common texts on numerical analysis.^ 

The accuracy of an answer produced by these algorithms is influenced by two 
factors: the performance of the subprogram (see Appendix B) and the accuracy 
of the argument. The effect of an argument error upon the accuracy of an answer 
depends solely upon the mathematical function involved and not upon the partic- 
ular coding used in the subprogram. 

A guide to the propagational effect of argument errors is provided because 
argument errors always influence the accuracy of answers whether the errors are 
accumulated prior to use of the subprogram or introduced by newly converted 
data. This guide (expressed as a simple formula where possible) is intended to 
assist users in assessing the effect of an argument error. 

The following symbols are used in this appendix to describe the effect of an 
argument error upon the accuracy of the answer: 

SYMBOL EXPLANATION 

q(x) The result given by the subprogram. 

i(x) The correct result. 

i(x) - g(x) 
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The relative error of the result 
fM given by the subprogram. 

8 The relative error of the argument. 

E I i(x) — g(x) I The absolute error of the result 

given by the subprogram. 

A The absolute error of the argument. 

The notation used for the continued fractions complies with the specifications 
set by the National Bureau of Standards.- 



iAny of the cominon numerical aiiahsis texts may be used as a reference. One such text is F. B. Hildebrand's 
Introduction In Sitincrwal Analysis (McGraw-Hill Book Company, Inc., N'ew York, X. Y.. 1956). Background 
information for alsxorithms that use continued fractions may be foinid in H. S. \\'airs Analytic Theory of 
Continued Fractions (I). VanXostrand Co., Inc., Princeton, X. J., 1948). 

-For more information, see Milton Abraniowitz and Irene A. Ste'.;nn (editors), Handl>ook of Mathematical 
Functions, Ajiplied Mathematics Serii's-So ( Xationa! Rin-ean of Standards, Washington, D.C.), 1965. 



xxxCLAES (CDABS) and xxxCSABS (CABS) Subprograms 

1. Write \x + iy\ = a + ib. 

2. lix = y = 0, then a = Oandb = 0. 

3. Let Vi = max ( |x|, |y| ), and 

V2 = min ( \x\, \y\). 



Then, a = v,'Jl+ (—Y, and ^ = 0. 



The algorithms for both complex absolute value subprograms are identical. 
Each subprogram uses the appropriate real square root subprogram (xxxlsqrt 

or 3CXXSSORT ) . 



xxxCLEXP (CDEXP) and xxxCSEXP (CEXP) Subprograms 
Algorithm 

The value of e-^+*J' is computed as e^ • cos(y) + i * e'^ * sm{y). The algorithms 
for both complex exponential subprograms are identical. Each subprogram uses 
the appropriate real exponential subprogram ( xxxlexp or xxxsexp ) and the appro- 
priate real sine/ cosine subprogram (xxxlscn or a:x%sscN). 

Effect of an Argument Error 

The effect of an argument error depends upon the accuracy of the individual 
parts of the argument, li e^+^y = R ' e'", then H = y and e{R) '- A (%). 

xxxCLLOG (CDLOG) and xxxCSLOG (CLOG) Subprograms 

Algorithm 

1. Write loge {x + iy) = a + ib. 

2. Then, a = loge \x + ij/jand b = the principle value of arctan - . 

The algorithms for both complex natural logarithm subprograms are identical. 
Each subprogram uses the appropriate complex absolute value subprogram 
(xxxcLAjBS or XXXCSABS ), the appropriate real natural logarithm subprogram 
(xxxLLOG or xxxslog), and the appropriate arctangent subprogram (xxxlatn2 

or XXXSATN2 ) . 

Effect of an Argument Error 

The effect of an argument error depends upon the accuracy of the individual 
parts of the argument, lix + iy = r ' e''' and loge {x + iy) = a + ib, then h = b 
andE(a) = 8{r). 

xxxCLSQT (CDSQRT) and xxxCSSQT (CSQRT) Subprograms 
Algorithm 



1. Write yx + iy — a -{■ ib. 

2. If X = y = 0, then a = Oandb = 0. 



3. If X > 0, then a = ^ — ■ — 

y 2 



and b — — 
2a 
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4. If X < 0, then a = -^ 
2a 



andZ? = (signy) '^— 



The algorithms for both complex square root subprograms are identical. Each 
subprogram uses the appropriate real square root subprogram (xxxlsqrt or 
xxxssqrt). 

Effect of an Argument Error 

The effect of an argument error depends upon the accuracy of the individual 
parts of the argument, li x + iy = r * &^ and yjx -\- iy = R • e^^, then 

e{R) -—S(r), and 6(H) ~8(/i). 



xxxCLSCN Subprogram (CDSIN and CDCOS) 
Algorithm 

1. If the sine is desired, then 

sin(x + iy) = sin(x) • cosh(j/) + i • cos(x) • sinh(j/). 
If the cosine is desired, then 

cos(x + iy) = cos(x) • cosh(i/) + i • sin(x) • sinh(j/). 

2. If X < 0, thensinh( — x) = — sinh(x). 

3. If x > 0=3465736, then sinh(x) = ^ 



2 
4. If ^ X ^ 0.3465736, then compute sinh (x) by use of the polynomial: 

cinn / v 1 

= flo + CliX^ + fl2X* + . . . + flsX^". 



X 

The coefficients are obtained by expanding the polynomial with respect to the 
Chebyshev polynomials over the range ^ x^ ^ 0.120113. The relative error 
of this approximation is less than 2~^^-^. 

5. The value of cosh(x) is computed as cosh(x) = sinhlx| + —^ . 

This computation uses the real exponential subprogram (xxxlexp) and the 
real sine/ cosine subprogram (xxxlscn). 

f ffecf of an Argument Error 

To understand the effect of an argument error upon the accuracy of the answer, 
the programmer must understand the effect of an argument error in the xxxlscn, 
XXXLEXP, and xxxlsnh subprograms. 



xxxCSSCN Subprogram (CSIN and CCOS) 

Algorithm 

1. If the sine is desired, then 

sin (x + iy) = sin(x) • cosh(t/) + i • cos(x) • sinh(y). 
If the cosine is desired, then 

cos(x + iy) = cos(x) • cosh(i/) + i* sin(x) • sinh(t/). 

2. If X < 0, then sinh( — r) = — sinh(x). 

1 

3. If X > 0.3465736, then sinh(x) = ^ . 



I 



20 



4. If ^ X ^ 0.3465736, then compute sinh (x) by use of the polynomial: 

sinh(.r) , ,, , , 
^—- ^ Qo + a^x- + a>x*. 

X 

The coefficients are obtained by expanding the polynomial with respect to the 
Chebyshev polynomials over the range ^ .^^ ^ 0.120113. The relative error 
of this approximation is less than 2~-''-^. 

5. The value of cosh(x) is computed as cosh(3c) = sinhja:! + — . 

This computation uses the real exponential subprogram (xxxsexp) and \he 
real sine/ cosine subprogram ( xxxsscn ) . 

Effect of an Argument Error 

To understand the effect of an argument error upon the accuracy of the answer, 
the programmer must understand the effect of an argument in the xxxsscn, xxxsexp, 
and xxxsscNH subprograms. 



xxxLASCN Subprogram (DARSIN and DARCOS) 

Algorithm 

1. If ^ X ^ — , then compute arccos(x) as: 

arccos ( x ) = -r — arcsin ( x ) . 
If ^ X ^ -^, then compute arcsin(x) by a polynomial of the form: 

Zi 

arcsin(x) = x + c^x^ + c-^x" + . . . + CriX-"^. 
The coefficients are obtained by expanding the function {{%) = ^—^, 

z = x^, with respect to the Chebyshev polynomials over the range, ^ x ^ -^. 
The relative error of this approximation is less than 2"^^-'^. 

2. If — < X ^ 1, then compute arcsin (x) as: 

2 ' '' 

If — < X ^ 1, then compute arccos (x) as: 

(1 1 — x\ 
A — -z — I . 

This case is now reduced to the first case because within these limits. 



0< llL£<± 
^=\ 2 = 2 ' 

This computation uses the real square root subprogram ( xxxlsqrt ) . 

3. If — 1 ^ X < 0, then arcsin(x) = — arcsin'xl 

andarccos(x) = tt — arccos |x|. 

This reduces these cases to one of the two positive cases. 
fffecf of on Argument Error 

E z===z . For small values of x, E -^ X Toward the limits ( ± 1 ) of the range 

VI - x^ 

a small A causes a substantial error in the answer. For the arcsine, e '--' 8 if the 

value of X is small. 
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IHCLATAN Subprogram (DATAN) 
Algorifhm 

1. Reduce the computation of arctan ( x ) to the case ^ x ^ 1 by using 
arctan(— x) = — arctan ( a: ) or 

arctan -j-nr = —^ arctan \x\ 

be 2 ' ' 



2. If necessary, reduce the computation further to the case Ixj ^ tan 15° by using 
arctan (x) = 30° + arctan 

V3'x - 1 



/ \/3'x-l \ 
\ x+ Y^ ) 



The value of 



^ tan 15° if the value of x is within the range, 



tan 15° < a: ^ 1. The value of ( V3 • x — 1) is computed as 
(V3 — l)x ^ ^-fxto avoid the loss of significant digits. 

3. For ]x\ ^ tan 15°, use a continued fraction of the form: 

arctan ( x ) _ ^ axX^ a^ a^ a^ 



X "~ (bi + x2) + (Z72 + x^) + (Z?3 + x") + (Z74 + x2) + ■ • • 

The relative error of this approximation is less than 2"^"^-^. The coefficients of 
this formula were derived by transforming the continued fraction: 

_J_ 3'4 16-25 

crofar. ( -y^ 3 25 * 7 7 ' 81 ' 11 

^^^_^^^^^ = 1 + 



a--o-(5%--o-G-^^^- 



4»3'9 
5 • 11 • 169 



1^13-17^^ / 



If 



where to has an approximate value oi ^^ ^ ~ „ {— x ~^ + 40) but the true 

64-27 

5j_289 • 19 

value of w is / yj^ 



\3-7-17 ) 



Effect of an Argument Error 

E-^ . For small values of x, e ^ 8, and as the value of x increases, the 

1 + x2 

efiFect of e upon 8 diminishes. 

XXXLATN2 Subprogram (DATAN and DATAN2) 

Algorithm 

IXi 1 7r 

-i- > 2^^, the answer = (sign Xi ) • -s- . 
X2 1 ^ 

Otherwise, if X2 > 0, the answer = arctan I — ) , and 

if X2 < 0, the answer = arctan ( — ) + ( sign Xi ) • tt. 
The rest of the computation is identical for either one or two arguments. 
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2. Reduce the computation of arctan(x) to the case 0^ ^c ^ 1, by using 
arctan(— x) = — arctan(x), or 






■~" arcLaniX 



2 
3. If necessary, reduce the computation further to the case jxj ^ tan 15° by using 

arctan(x) = 30° + arctan z=^. 

x+ y/3 

V3 ■ X - 1 

X + V 3 

tan 15° < X ^ 1. The value of ( y/3 • x — 1) is computed as ( ^/3 — 1) x — 1 



The value of 



^ tan 15° if the value x is within the range, 



to avoid the loss of significant digits. 
4. For |x| ^ tan 15°, use a continued fraction of the form: 

arctan(x)_^ a^x^ ^2 as a4: 



The relative error of this approximation is less than 2"^'^-^. The coefficients of 
this formula were derived by transforming the continued fraction: 

j^l 3*4 16*25 

arctan(x) ^ ^ _^ 3 25 * 7 7 * 81 • 11 



(I--)- (^--')- (s^3^-0 



4*3*9 
5 • 11 • 169 



-2 I _ 



\13*17 ^ ) 



w 



2 



where to has an approximate value of ■= — — — r-:; — t;= ( —x~'^ + 40) but the 

5 * 11 * 13 * 17 ^ 

64*27 

1 r . 5*289*19 

true value of w is 



\3*7*17 )^ 

Effect of an Argument Error 

E '— -p— — r . For small values of x, e '-' S, and as the value of x increases, the effect 
1 + x-^ 

of e upon 8 diminishes. 

xxxLERF Subprogram (DERF and DERFC) 
Algorithm 

1. If ^ X < 1, then compute the error function by the following approximation: 

erf(x) ^ X (flo + dix^ + azX^ + ., . + a-nX^^). 

Grf i Oil 
The coefficients were obtained by expanding the function {{z) — ^-^, 

z — x^, with respect to the Chebyshev polynomials over the range, ^ x < 1. 

The relative error of this approximation is less than 1.07 * 2~^^. The value of 

the complemented error function is computed as erfc(x) = 1 — erf(x) and is 

1 
greater than r- . 
Id 

2. If 1^ X < 2.0400009, then compute the complemented error function by the 
following approximation: 

erfc (x) =bo + hiZ + h-iZ^ + . . . + h^^z^^ 
where z = x — Tq and To = 1.999999if;. The coefficients were obtained by ex- 
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panding the function i{z) = erfc(z + To) with respect to the Chebyshev poly- 
nomials over the range — 1 ^ x ^ 0.04. The absolute error of this approxima- 
tion is less than 1.5 • 2^*^^ The limits of this range and the base value for To 
were used to minimize the hexadecimal truncation error. The value of the com- 
plemented error function within this range is greater than -^^ . The value of 
the error function is computed as erf (x) = 1 — erfc(x). 

3. If 2.0400009 ^ X ^ 13.306, then compute the complemented error function by 

the following approximation: 

^ . , (co + CiX-2 + C2X-4 + . . . + C2oX-^o) e-^^ 

ertc ( X ) = -^^ '■ 

^ ' X 

The coefficients were obtained by expanding the function 
i{z) = erfc (x) • x • e^% z = x-^, with respect to the Chebyshev polynomials 
over the range 2.04"- > z ^ 13.306"^. The relative error of this approximation 
ranges from 2~^^ at 2.04 to 2-^^ at 13.306. This computation uses the real ex- 
ponential subprogram ( xxxlexp ) . 

If X ^ 6.092, then the error function is computed as erf (x) = 1 — erfc(x). 
If X > 6.092, then the error function is ^ 1. 

4. If 13.306 < X, then the error function is ^ 1, and the complemented error 
function is ^ 0. 

5. If X < 0, then reduce to a case involving a positive argument by the use of 
the following formulas: 

erf(— x) = — erf(x) and erfc(— x) = 2 — erfc(x). 



E ^ e~^^ * A. For the error function, as the magnitude of the argument exceeds 
1, the effect of an argument error upon the final accuracy diminishes rapidly. For 
small values of x, e '^ S. For the complemented error function, if the value of 

X is greater than 1, erfc(x) '^ - . Therefore, c ^ 2x^ • S. If the value of x 

is negative or less than 1, then e '^ e~'^^ • A. 

xxxLEXP Subprogram (DEXP) 

Algorlihm 

1. If X < — 180. 2183, then is given as the answer. 

2. Divide x by loge2 and write 

w = = (4a — h d) 

loge2 16 

where a, b, and c are integers, ^ fo ^ 3, ^ c ^ 15, and d is within the range 

< d < -i_ . Then, e^ = 2^ = 16« • 2-* • 2-"^^^ • 2-". 
~ 16 

3. Compute 2~^ by using the Chebyshev interpolation of degree 6 over the range, 
= <^ < ■ The maximum relative error of this computation is 2~^^. 

4. If c > 0, then multiply 2-^ by 2-«/i6. (The 15 values of 2-'^^^ for 1 ^ c ^ 15 
are included in the subprogram. ) 

5. If Z? > 0, then halve the result b times. 

6. Finally, add the hexadecimal exponent a to the characteristic of the answer. 

Effect of an Argument Error 

E '-' A. If the magnitude of x is large, even the round-off error of the argument 
causes a substantial relative error in the answer because A = 8 • x. 



I 
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xxxLGAMA Subprogram (DGAMMA and DLGAMA) 

Algorithm 

1. If < X ^ 2~-^'-, then compute log-gamma as logpr(x) = — logc(x). This 
computation uses the real logarithm subprogram ( xxxllog ) . 

2. If 2" -•"'•- < 1' < 8, then compute log-gamma by taking the natural logarithm 
of the value obtained for gamma. The computation of gamma depends upon 
the range into which the argument falls. 

3. If 2~-^'~ <C X <C I, then use r(x) = — ^^ — - to reduce to the next case. 

X 

4. If 1 ^ X ^ 2, then compute gamma by the following approximation: 

r{x) ^ Go + fllZ + (I2Z' + . . . + ^22^^^ 

where z = x — 1.5. The coefficients were obtained by expanding the function 
f{z) = T{x) with respect to the Chebyshev polynomial for jz| ^ 0.5. The 
absolute error of this approximation is less than 1.5 • 2~^^. 

5. If 2 < X < 8, then use T{x) = (x — I) Y{x — 1) to reduce to the preced- 
ing case. 

6. If 8 ^ X, then compute log-gamma by the use of Stirling's formula: 

logj(x) =x(l0g,(x) - 1) - ^loge(x) + -|-l0ge(27r) + G(x). 

The modifier term G(x) is computed as 

G(x) = Z?ax-i + b.>x-'^ + bsx-'' + b^x-' + b-^x-\ 

G(x) 
The coefficients were obtained by expanding the function f{z) = — ^-^ , 

z = x~-, with respect to the Chebyshev polynomials over the range 
< z < 8~-. The absolute error of the approximation for G(x) is less than 
X • 2"^''*''. Because, in this range, x < loger(x), the contribution of this error 
to the relative error of the value for log-gamma is less than 2~^^. This com- 
putation uses the real logarithm subprogram ( xxxllog ) . 

For gamma, compute r(x) = e'-', where y is the value obtained for log- 
gamma. This computation uses the real exponential subprogram (xxxlexp). 

Effect of an Argument Error 

£ -^ "/'(^) * ^ for gamma, and E r^ i^i^) * ^ for log-gamma, where i/' is the 
digamma function. 

If -^ < X < 3, then —2 < i}/{x) < 1. Therefore, E ^ A for log-gamma. 

However, because x = 1 and x = 2 are zeros of the log-gamma function, even 
a small 8 can cause a substantial e in this range. 

If the value of x is large, then i/'(x) ~ loge(x). Therefore, for gamma, 
e ^ 8 • X • loge(x). In this case, even the round-off error of the argument con- 
tributes greatly to the relative error of the answer. For log-gamma with large 
values of X, € ^ S. 

xxxLLOG Subprogram (DLOG and DLOG10) 

Algorithm 

1. Write X = 16'' • 2~^ • m. where p is the exponent, q is an intef^er, ^ 't ^ 3, 

and m is within the range, — < m < 1. 

2 ~ 

2. Define two constants, a and b (where a = base point and 2~^ = a) as follows: 

If -v" — m < ^=, then a = -r- and b = 1. 

-^ - \/2 2 

If — = < m < 1, then a = I and b = 0. 

\'2 ~ 
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tn — n 1 -I- ? 

3. Write z = , . Then, m = a • , and |z| < 0.1716. 

m + a 1 — z 

1 + z /l + zX 

4. Now, X = 24P-9-6 •-——-, and log^x = {4p - q - b) loge2 + loge I , _ j . 

5. Finally, loge ( :j I is computed by using the Chebyshev interpolation of 

degree 7 in z^ over the range, O^z^^O. 02944. The maximum relative error 
of this approximation is 2"~°^-^. 

6. If the common logarithm is desired, then logioX = logioe • loge^c. 

Effecf of an Argument Error 

E '— ' 8. Therefore, if the value of the argument is close to 1, the relative error can 
be very large because the value of the function is very small. 

xxxLSCN Subprogram (DSIN and DCOS) 
Algorithm 

1. Divide |x| by-^- and separate the quotient (z) into its integer part (q) and 

4 
its fraction part (r). Then, z = x • — = q + r, where q is an integer and r 

is within the range, ^ r < 1. 

2. If the cosine is desired, add 2 to (7. If the sine is desired and if x is negative, 
add 4 to q. This adjustment of q reduces the general case to the computation of 
sin (x) for x ^ 0, because 

cos ( ± X ) = sin I |xl + ^ J , and 

sin ( — x) = sin ( |x| + tt). H 

3. Let qo^q mod 8. fl 

Then, for go = 0, sin ( x ) = sin I -j- • r I 

qo = 1, sin (x) = cos ( "|"(1 ~ ^) ) 

go = 2, sin ( X ) = cos ( "J" * ^ I 

go = 3,sin(x) = sin (^(1 - r)\ 

go = 4, sin (x) = —sin ( "7"* ^ I 

go = 5, sin (x) = —cos I -j~{^ ""'')) 

go = 6, sin (x) = —cos I "7~* ^ I 

go = 7, sin (x) = -sin (^(1 " *") I 
These formulas reduce each case to the computation of either sin 1 -j-* fi I or 
cos I -Y' ^1 I ; where ri is either r or (1 —r), and is within the range, O^fi ^ 1. 

4. Finally, cither sin j -^* r^ \ or cos J -j-* fi j is computed, using the Chebyshev 

interpolation of degree 6 in fi^ for the sine, and of degree 7 in rr for the cosine. 
The maximum relative error of the sine polynomial is 2""^^ and that of the 
cosine polynomial is 2~^^-^. 
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Effect of an Argument Error 

\'alue diminishes periodically, no consistent relative error control can be main- 
tained outside of the principal range, ^ ^ -'^' ^ + ~5~ ■ 



xxxLSCNH Subprogram (DSINH and DCOSH) 

Algorithm 

1. If .V < 0.3465736, then compute sinh(x) as: 

„:„U / -, \ ^, .. I „ ..-A 1 „ ..r> I „ ..7 I „ „>) 1 „ ,.11 

SHJll ■ A I == A ~r L-;A ~r L'-jA "T C.vA " HT C'4A T~ '--5A '■ 

sinh ( Xl 

The coefRcients are obtained by expanding the function f{z) = ^-—^ , 

z = X-, with respect to the Chebyshev polynomials over the range, 

^ ;: < 0.12011326. The relative error of this approximation is less than 2""'^'". 

2. If either ^r ^ 0.3465736 or the cosh(x) is desired, obtain if = e''L Then, 
cosh(.t) = , and sinh(x) = (sign x) • . The real expo- 
nential subprogram (xxxlexp) is used to compute the value of tc. 

Effect of an Argument Error 

For the hyperbolic sine, E ~ A • cosh(x) and e ~ A • coth(x). 
For the hyperbolic cosine, E ~ A • sinh(x) and e ~ A • tanh(x). 

Specifically, for the cosine, E ~ A over the entire range; for the sine, e ^ 8 
for the small values of x. 



xxxLSQRT Subprogram (DSQRT) 

Algorithm 

1. If X = 0, then the answer is 0. 

2. Write x = 16-^~^ • m, where 2p — q is the exponent and q equals either or 1; 
m is the mantissa and is within the range, — ^ ?7i ^ 1. 

3. Then, V^ = 16^ ' S"-*^ • \hn . 

4. For the first approximation of \'x , compute the following: 

i/o = 2--^-16^-('-|--f-^.mV 

The maximum relative error of this approximation is — . 

5. Apply the Newton-Raphson iteration 

four times to ?/o (the first two times in the short form and the last two times in 
the long form). The final step is performed as 

^^ = ^^^ + t(-^ - ^0 

to minimize the computational truncation error. The maximum relative error 
of the final result is theoretically 2~''^'''^\ 

Effect of an Argument Error 

1 



8 
9 
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xxxLTANH Subprogram (DTANH) 
Algorithm 

1. If \x\ < 0.54931, then use the following fractional approximation: 
tanh (x) _ ^ _ flix - + o.^^c* + fl:)^;^ + x^ 



X hi) + hix- + Z?^^'' + ^3^'^ + x^ 

where: 

ai = 676440.765 Z?o = 2029322.295 

fl2 = 45092.124 hi = 947005.29 

as = 594.459 h. = 52028.55 

hs = 630.476 

The maximum relative error t)f this approximation is 2~*'^-'"'. The formula was 
obtained by transforming the continued fraction 
tanh (x) _ 1 , x^ x^ x^ 



x 3+5+""*15 + u; 

where w has an approximate value of 0.017, but the true value of to is 
x^ x^ 
ITT 19+ • * ■ 

2 

2. If 0.54931 ^ X < 20.101, then use the identity tanh (x) = 1 iT^TY" 

computation uses the double precision exponential subprogram (xxxlexp). 

3. If X ^ 20.101, then tanh (x) s 1. 

4. If x^ — 0.54931, then use the identity tanh (x) = — tanh ( — x). 

Effecf of an Argument Error 

2 A 
E -^ ( 1 — tanh^ x) A, and e '^ .,,-.. . For small values of x , £ ^ 8. As the 

smh (2x) 

value of X increases, the effect of S upon e diminishes. 



xxxLTNCT Subprogram (DTAN and DCOTAN) 
>\/gorif/iin 

1. Divide |xj by ^ and separate the result into the integer part (q) and the 

fraction part (r). Then, |x] = -r- {q + r). 

2. Obtain the reduced argument {w) as follows: 

if q is even, then w = r. 
if q is odd, then «; = 1 — r. 
The range of the reduced argument is ^ u) ^ 1. 

3. Let qo = q mod 4. 

Then, for qo = 0, tan |x[ = tan i~7^'^ I and cot |x| = cot I -j- • if j 

<7o = 1, tan |x| = cot I ~7" * «^ I and cot [xj = tan I "X * ^^ ) 
qo = 2, tan [x| = —cot ( ^ * "^ I and cot [x| = —tan ( ~7" * ^^" J 
qo = 3, tan |x[ = —tan ( ~7~ * w; j and cot jx] = —cot I ~:r * ^^^ ) 

4. The values of tan ( -7- • to\ and cot I -j" * ii^\ are computed as the ratio 
of two polynomials. 



I 
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foj-j f .. — • -»/' I --^^ i L r^Tr\r\ r\r\¥ f -11 • ^j* 1 ^—^^ t \ / 



(t* ^0 "'■^^^^^''"'^'°'(^' "") 



where P{ic-) is of degree 3 and Q{tv-) is of degree 4 in lo^. The coefficients 
of P and Q are obtained by economizing the continued fraction 

tan(;::) £^ _^ ^ 

~ X fy ^ p^ . . . 

~ 3— o— i — 

in the following way. 

Write: ^-^^^^ ^ i - -f_ Ji_ -^ 



3- 5- 7- 9- (11 + rfi)- (13 + ^2)- (15 + ^3) 
and determine the values for d^, d->, and d-^ so that the right-hand expression 

relative error of this formula over the range < z < -^ is 3.4 • 10"^-'. 

o — — 4 

Change the variable from z io w = -^ • z and rewrite the formula to obtain 

4 

P{\c-) and Qlyxc-). 
5. If X < 0, then tan(x) = —tan \x\,, and cot(x-) = —cot \x\. 

Effect of an Argument Error 

A , 2 . , - 

E ~ — —. — : , and e ~ . ,. . for tan (' x ) . Therefore, near the singularities 
cos-(x) sm(2x) ^ ^ ^ 

ofx=llc+-r-l7r, where k is an integer, no error control can be maintained. 
This is also true for cotan(.i:) for values of x near kv, where k is an integer. 

xxxSASCN Subprogram (ARSIN and ARCOS) 
Algorithm 

1. If ^ X ^ — , then compute arccos(x) as: 

arccos(x) =— — arcsin(x). 

If ^ X ^ — , then compute arcsin(x) by a polynomial of the form: 

arcsin(x) ^ x + CyX^ + c-vx^ + Csx'^ + c^x^ + CsX^'. 

3.rcsin ( x ) 
The coefficients are obtained by expanding the function f(z) = ^^ — -, 

z — X-, with respect to the Chebyshev polynomials over the range ^ z ^-r-. 
The relative error of this approximation is less than 2~-^ •''. 

2. If — < X ^ 1, then compute arcsin(x) as: 

arcsin(x) = -z — arccos(x). 
If -T- < X ^ 1, then compute arccos(x) as: 

arccos(x) = 2 • arcsin l A — - — j. 
This case is now reduced to the first case because within these limits, 



11 - X 1 



I 



2 = 2 ■ 

This computation uses the real square root subprogram ( xxxssqrt ) . 

3. If — 1 ^ X < 0, then arcsin(x) = — arcsin !x], and arccos(x) = tt — arccos |x|. 
This reduces these cases to one of the two positive cases. 
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Effect of an Argument Error 

E '■^ —= . For small values of x, E -^ A. Toward the limits ( ± 1 ) of the 

VI — de- 
range, a small A causes a substantial error in the answer. 



IHCSATAN Subprogram (ATAN) 
Algorithm 

1. Reduce the computation of arctan ( x ) to the case ^ x ^ 1, by using 

arctan ( —x)-= —arctan (x), or 

arctan I TT I ~ ^1 — arctan |x|. 

2. If necessary, reduce the computation further to the case |x[ ^ tan 15° by using 

arctan (x) = 30° + arctan / V3 ' x - 1 \ ^ 

\ X + V3 / 



The value of 



\/3*x - 1 



^ tan 15° if the value of x is within the range. 



X + V3 

tan 15° < X ^ 1. The value of (V3 • x — 1) is computed as 
( V3 — l)x — 1 + xto avoid the loss of significant digits. 
3. For [x[ ^ tan 15°, use the approximation formula: 

5^^^2!LM^ 0.60S10579 - 0.0516O454x' + J'f"^,, 

X x^ + 1.4087812 

±iiib luiiiiuici iicti a. icictuvc ciiui icaa iiictii .ij ttiiu can uc uuLctiiicLi uy 

transforming the continued fraction 

arctan (^c) _ , _ x^ 5 



3 + / 5 

' ' — w 



— — x~- H — ==- 1 10~*, but the true 



I 



4'5 
7.7.9 
value of w is 



(7^-i 



+ 



The original continued fraction can be obtained by transforming the Taylor 
series into continued fraction form. 



Effect of an Argument Error 

A. 
1 + x 
of 8 upon € diminishes 



^ '^ T _L_ „2 • For small values of x, € ^ 8; as the value of x increases, the effect 



XXXSATN2 Subprogram (ATAN and ATAN2) 

Algorithm 

1. For arctan(X], X2), if either X2 = or — ^ > 2^^, the answer = (sign Xi ) • -r-. 

I X2 1 '^ 

Otherwise, if Xj > 0, the answer = arctan 1 — 1 , and 

\xj 

if X2 < 0, the answer = arctan ( -^ I + ( sign -^i ) ^r. 

The rest of the computation is identical for either one or two arguments. 
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2. Reduce the computation of arctan(x) to the case ^ x ^ 1, by using 
arctan(— x) = — arctan(a:) or 



arctan l -r-r I =-:;^ — arctan Ixj. 



(w) 2 
luce the compu 

arctan(x) = 30° + arctan- 



The value of 



^ tan 15° if the value of x is within the range, 



3. If necessary, reduce the computation further to the case jx| ^ tan 15° by using 

V3'x- 1 

x+ V3 ■ 
V'3 • ac — 1 

x+ V3 
tan 15° < X '^ 1. The value of (V3 • ^ ~ 1) is computed as 
( V3 — l)x— 1+xto avoid the loss of significant digits. 

4. For [x[^ tan 15°, use the approximation formula: 
arctan(x) _ o_o60310579 - 0.0516045x2 + 0.55913709^ 



^ ...^.v...^^ ■ x2 + 1.4087812 

This formula has a relative error less than 2"^'^-^ and can be obtained by trans- 
forming the continued fraction 



arctan(x) _ -i _ ^^ 5 

~ 3 + / 5 

■ " w 



— 77 ^~^ + ~77~ ) ^^~* ^^* *^^ ^^^ 



4'5 

7.7.9 
value of IV is 



(^-i 



+ 



The original continued fraction can be obtained by transforming the Taylor 
series into continued fraction form. 

f fFecf of an Argument Error 

E r^ - — ; — 5 . For small values of x, e --^ 8; as the value of x increases, the effect 
1 + x^ 

Oi o Upon e u-imimsucS. 



xxxSERF Subprogram (ERF and ERFC) 
Algorithm 

1. If ^ X ^ 1.317, then compute the error function by the following approxi- 
mation: 

erf(x) = x(ao + ^i^^ + (^2^^ + . . . + a^x^^). 

Grr ( Xt 
The coefficients were obtained by expanding the function f ( z ) = — ^— ^ ,z = x^, 

with respect to the Chebyshev polynomials over the range ^ x ^ 1.317. The 
relative error of this approximation is less than 2-24_ xhe value of the comple- 
mented error function is computed as erfc(x) = 1 — erf(x) and is greater 

than-^. 

2. If 1.317 < X ^ 2.0400009, then compute the complemented error function by 
the following approximation : 

erfc(x) ^bo + hxZ + Z?2Z" + . . . + h-iZ^ 
where z = x — To and To = 2.0400009. The coefficients were obtained by 
expanding the function f(z) = erfc(x + Tq) with respect to the Chebyshev 
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polynomials over the range ( 1.317 — To) < z ^ 0. The absolute error of this 
approximation is less than 1.3 • 2~'^*'. The value of the complemented error 

function within the range 1.317 < x ^ To is greater than -z;^ . The value of the 
error function is computed as erf (x) = 1 — erfc(x). 

3. If To < ic ^ 13.306, then compute the complemented error function by the fol- 
lowing approximation: 

r , . (Co + CiX-^ + C^X--^ + . . . + C6X-12\ q-x"^ 

ertc(x) ^^- 

X 

The coefficients were obtained by expanding the function 
f(z) = erfc {x) • X* ^\z — x~^, with respect to the Chebyshev polynomials 
over the range Tq"^ > ^ ^ 13.306^^. The relative error of this approximation 
is less than 1.2 • 2"^^. This computation uses the real exponential subprogram 

( XXXSEXP ) . 

If X ^ 3.9192, then the error function is computed as erf (x) = 1 — erfc(x). 
If X > 3.9192, then the error function is ^ 1. 

4. If 13.306 < X, then the error function is ^ 1, and the complemented error func- 
tion is ^ 0. 

5. If X < 0, then reduce to a case involving a positive argument by the use of the 
following formulas: 

erf(— x) = — erf(x) and 
erfc(— x) = 2 — erfc(x). 

Effect of an Argument Error 

E '^ e"'^^ • A. For the error function, as the magnitude of the argument exceeds 1, 
the effect of an argument error upon the final accuracy diminishes rapidly. For 
small values of x, e -~' 8. For the complemented error function, if the value of x is 

greater than 1, erfc( x ) '-' -r— . Therefore, e ^ 2 x^ • 8. If the value of x is negative 
or less than 1, then e ^ e~^'^ • A. 

xxxSEXP Subprogram (EXP) 
A/gor/ffim 

1. If X < — 180.218, then is given as the answer. 

2. If [x] < 2~2^, then 1 is given as the answer. 

3. Otherwise, divide x by loge2 and write 

where a and h are integers, ^ & ^ 3 and ^ d < 1. 
Then, e^ = 2«' = (16« • 2"* •2-'^). 

4. Compute 2"*^ by the following fractional approximation: 



2-d^ 



0.034657359 d^ + d + 9.9545958 - ^^'^■^'^227 



d^ + 87.417497 



This formula can be obtained by transforming the Gaussian-type continued 
fraction 

^ _ , z z z z z z z z_ 

^ ~ 1+2- 3+ 2- 5+ 2-7+2 

The maximum relative error of this approximation is 2~^^. 

5. Multiply 2-"^ by 2-^ 

6. Finally, add the hexadecimal exponent a to the characteristic of the answer. 
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Effect of an Argument Error 

e -^ A. If the magnitude of x is large, even the round-off error of the argument 
causes a substantial relative error in the answer because A = S • x. 

xxxSGAMA Subprogram (GAMMA and ALGAMA) 
Algorithm 

1. If < X ^ 2--^~, then compute log-gamma as loge T (x) = —lege (x). This 
computation uses the real logarithm subprogram ( xxxslog ) . 

2. If 2"^^^ < X < 8, then compute log-gamma by taking the natural logarithm of 
the value obtained for gamma. The computation of gamma depends upon the 
range into which the argument fails. 

3. If 2~2^^ < X < 1, then use F (x) = — ^^ to reduce to the next case. 

4. If 1 ^ X ^ 2, then compute gamma by the following approximation: 

r (x) = flo + «l2 + a2Z^ + . . . + flgZ^ 

where % = x — 1.5. The coefficients were obtained by expanding the function 
f(z) = r (x) with respect to the Chebyshev polynomials for \z\ ^ 0.5. The 
absolute error of this approximation is less than 1.5 • 2~^^. 

5. If 2 < X < 8, then use r(x) = {x -- l)r(x— l)to reduce step by step to 
the preceding case. 

6. If 8 ^ X, then compute log-gamma by the use of Stirling's formula: 

I0ger(x)-X(l0ge(x) - 1 ) " y loge ( X ) + |- log,. ( 2x) + G( X ) . 

The modifier term G(x) is computed as 

G(x) = foix-i -1- Z?2X-2. 
The absolute error of the approximation for G(x) is 1.4 • 2~^-^ This computa- 
tion uses the real logarithm subprogram ( xxxslog ) . 

For gamma, compute F (x) = e^', where y is the value obtained for log- 
gamma. This computation uses the real exponential subprogram ( xxxsexp ) . 

Effect of an Argument Error 

£ '^ i/^ (x) • A for gamma, and E ^ i// (x) • A for log-gamma, where ij/ is the 
digamma function. 

If — < X < 3, then — 2 < i/^ (x) < 1. Therefore, E '-' A for log-gamma. How- 
ever, because x = 1 and x = 2 are zeros of the log-gamma function, even a small 
8 can cause a substantial e in this range. 

If the value of x is large, then \fr (r) ~ log^ (x). Therefore, for gamma, 
e ~ 8 X • loge (^)- In this case, even the round-off error of the argument con- 
tributes greatly to the relative error of the answer. For log-gamma with large 
values of X, c '^ 8. 

xxxSLOG Subprogram (ALOG and ALOG10) 
Algorithm 

1. Write X = 16^ • m, where p is an integer and m is within the range, -r^ ^ m < 1. 

2. Define two constants, a and b, where a = base point and '2~'' = a, as follows: 

If -T^ ^ m < -r-, then a = -rr- and b = 4. 

ID O ID 

If -3- ^ m < -r-, then a = —r- and b = 2. 
o z 4 

If -r- ^ m < 1, then a =1 and b = 0. 
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3. Write z = -^^^— ; . Then, m = a* —^ , and jz] ^ -75- . 

m + a 1 — z o 

4. Now, X = 2*P-^' ^ _ ^ , and logeX = (4p - b) loge2 + log^ I , _ ^ I . 



1 - 



5. Finally, loge ( _ I is evaluated using the Chebyshev interpolation of degree 

4 in z~ over the range, ^ z^ ^ -rr- . The maximum relative error of this ap- 
proximation is 2""^^-^. 

6. If the common logarithm is desired, then logioX = logioe • logeX. 

Effect of an Argument Error 

E '-' 8. Specifically, if 8 is the round-off error of the argument, e.g., 8 '-' 6 • 10~^, 
then E -^ 6 • 10~^. Therefore, if the argument is close to 1, the relative error can 
be very large because the value of the function is very small. 



xxxSSCN Subprogram (SIN and COS) 
Algorithm 

4 

1. Define z = — • x and separate z into its integer part (q) and its fraction part 

(r).Thenz = q + r, and Ix| = I ^ • <? I + ( ~7"* ^ ) ■ 

2. If the cosine is desired, add 2 to q. li the sine is desired and if x is negative, 
add 4toq. This adjustment of q reduces the general case to the computation 
of sin (x) for X ^ because 

cos ( ± X ) = sin I -|- + X I , and 

sin ( — x) = sin (tt + x). 

3. Let qo ^ q mod 8. 

Then, for qo = 0, sin (x) = sin I -^* r j 

qo = l,sin(x) = cos 1-^(1 - r) \ 

(7o = 2, sin (x) = cos I -7- • r I 

qo = 3, sin(x) = sin (^ (1 - r) \ 

qo = 4, sin (x) = -sin ( "^' ^1 
qo = 5, sin (x) = —cos 1 ~r{^ ~ ^) ) 
qo = 6, sin (x) = —cos I -j-* r I 
qo = 7,sin(x) = -sin ("f-(l ~ ^) ) 
These formulas reduce each case to the computation of either sin I — • r^ I or 

cos ( X* ^1 I where n is either f or (1 — r) and is within the range, 
< fi < 1. 
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4. Finally, the computation for either the sine or the cosine is performed, using 
the Chebyshev interpolation of degree 3 in r-^. The maximum relative error 
of the sine polynomial is 2~^^-^ and that of the cosine polynomial is 2~2^-^. 

Effect of an Argvment Error 

E -^ A. As the value of x increases, A increases. Because the function value dimin- 
ishes periodically, no consistent relative error control can be maintained outside 

the principal range, — «-^ x ^ + -^ . 

xxxSSCNK Subprogram (SINH and COSH) 
Algorithm 

1. If |x[ < 0.3465736, then compute sinh(x) as: 

sinh(x) ^ X + 0.16666505x3 + 0.00836915x5. 

The coeflBcients were obtained by expanding the function i{z) = ^, 

z = x^, with respect to the Chebyshev polynomials over the range 
< z < 0.12011326. The relative error of this approximation is less than 

2 — 26.5 

2. If either |x| ^ 0.3465736 or the cosh(x) is desired, obtain w = gl^L Then, 
cosh(x) = r , and sinh(x) = (sign x). ^ . The real exponen- 
tial subprogram ( xxxsexp ) is used to compute the value of to. 

Effect of an Argument Error 

For the hyperbolic sine, E '^ A • cosh(x) and e '^ A* coth (x). 
For the hyperbolic cosine, E --w A • sinh (x) and e '~ S • tanh(x). 

Specifically, for the cosine, c /^ A over the entire range; for the sine, e ^ 8 for 
small values of x. 



xxxSSQRT Subprogram (SQRT) 
Algorithm 

1. If X = 0, then the answer is 0. 

2. Write x = 16^^ • m, where p is an integer and m is within the range, 

2k=^<^- 

3. Then, V^ = 16^ • \/m., where p is the exponent of the answer and m is the 
mantissa of the answer. 

4. For the first approximation of 's/m, take hyperbolic approximations of the form 

b 

a -\ — — where the values of a, b, and c depend upon the value of m as follows: 

a. If -j^ ^ m < 1, then a = 1.80713 
b= -1.57727 
c = 0.954182 

These values minimize the maximum relative error (eo) over the range, 
while making an exact fit at m = 1. The exact fit at m = 1 minimizes the 
computational loss of the last hexadecimal digit for the values of m slightly 
less than 1. The relative error of this approximation is less than 2~^-^^. 
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b. If -^ ^ m <^ , then a = 0.428795 
^^^ ^^ b= -0.0214398 

c = 0.0548470 

These values minimize m^^^ • eo over this range of m where co denotes the 
relative error of this approximation, to is less than 2~^-^ • m"^/®. 

5. Multiply the result by 16*" to obtain the first approximation (i/o) of the answer. 

6. To obtain the final answer, the Newton-Raphson iteration 



H'-^i) 



\ f/n / 

must be applied twice to t/o- For — ^ m < 1, the final relative error is 

theoretically less than 2-2^'^; for -^p^ ^ m <Cjw, the final absolute error is 
theoretically less than 2-^9 • 16^. 

Effect of an Argument Error 

1 



xxxSTANH Subprogram (TANH) 
Algorithm 

1 Tf UI <r 0-12 4-Kor> fonV. f^\^ 

JL» XJL UV I — -I ^J ^ l.AJ.^/J.X CCAXXJLJL \ -^ y 



2. If 2~^2 < |x| < 0.54931, use the following fractional approximation: 
tanh (x) ^ t" + 3 5.1535 



2 . ^sic^o _L 105.4605 
x^ + 45.1842 +■ 



x- 

This approximation has a relative error less than 2~^^. The formula can be 
obtained by transforming the continued fraction 
tanh (x) _ 1 , £^ _^ x^ 



X 3+ 5+ 7+u; 

where w has an approximate value of 0.0307, but the true value of to is 

9+ llT"'* 

3. If 0.54931 ^ X < 9.011, then use the identity tanh (x) = 1 - ^^^^ y The 
computation for this case uses the real exponential subprogram (xxxsexp). 

4. If X ^ 9.011, then tanh (x) ^ 1. 

5. If x^ —0.54931, then use the identity tanh (x) = — tanh(— x). 

Effect of an Argument Error 

2a 
E ^ ( 1 — tanh^ x ) A, and e r^ — .,,-.. . For small values of x, e '--' 8, and as the 

smh (2x) 

value of X increases, the effect of 8 upon e diminishes. 

xxxSTNCT Subprogram (TAN and COTAN) 
Algorithm 

1. Divide |x| by "2- and separate the result into the integer part {q) and the 

fraction part (r). Then, [x| = -7^ (Q + ^)- 
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2. Obtain the reduced argument (it;) as follows: 

if q is even, then vj = r, 
if q is odd, then w = I — r. 
The range of the reduced argument is ^ it; ^ 1. 

3. Let qo^q mod 4. 

Then, for qo = 0, tan \x\ = tan ( "7" * *^ ) ^^^ ^^^ !^1 = cot j -^ • if I 
qo = 1, tan |a:| = cot I -j- • if j and cot \x\ = tan ( ~j" * «^ I 

ri^ zz: V- ran w\ ^ — r-^nr a . •-!/■; a o-nn *^rir i-v-i ;^: — t-on a • -ff- i 

q'o = 3, tan jx[ = —tan ( "j" ' w^ I and cot \x\ = —cot I -j- • w j 

4. The values of tan I ~7~ * w; J and cot I ^ • w\ are computed as the ratio of 
two polynomials. 

tan I -J- • 10 I s ^^ -, and cot I -r * w? I = ^ — 

where P{w^) = 212.58037 - 12.55991210^ 

Q(«;2) = 270.665736 - 71.645273w;2 + m;" 
This approximation is obtained by economizing the continued fraction 
tan(z) _ z^ z^ z^ 

z ^ 3-5-7-'"" 

in the following way: 

Write: ^^^^^1- ^' 



z - (3 + di)- (5 + ^2)- (7 + ^3) 
and determine values for di, 3,2, and d^ so that the right-hand expression 
gives the exact answers for z^ = 0.19, 0.432, and 0.594. Then the maximum 

relative error of this formula over the range ^ z ^-j-is 1.74 • 10~^. 

4 
Change the variable from z to to = — • z and rewrite the formula to ob- 

tain P(w^) and Odj^^) 
5. If X < 0, then tan \x\, and cot(x) = —cot |x|. 

Effect of an Argument Error 

A 2 

E '^ — 57~~^ , and e ~ . ,^ . for tan(x). Therefore, near the singularities 
cos^(x) sm(2x) ^ ^ * 

X = ik +-^J7r, where k is an integer, no error control can be maintained. 
This is also true for cotan(x) for x near k-n-, where k is an integer. 
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Appendix B. Performance Statistics 



Appendix B contains accuracy and timing statistics for 
the explicitly called mathematical subprograms. These 
statistics are presented in Table 12 and are arranged 
in alphabetical order, according to the entry names. 
The following column headings are used in Table 12: 

Entry Narhe: This column gives the entry name that 
must be used to call the subprogram. 

Argument Range: This column gives the argument 
range used to obtain the accuracy figures. For each 
function, accuracy figures are given for one or more 
representative segments within the valid argument 
range. In each case, the figures given are the most 
meaningful to the function and range under consid- 
eration. 

The maximum relative error and standard deviation 
of the relative error are generally useful and revealing 
statistics; however, they are useless for the range of a 
function where its value becomes 0, because the slight- 
est error in the argument can cause an unpredictable 
fluctuation in the magnitude of the answer. When a 
small argument error would have this effect, the maxi- 
mum absolute error and standard deviation of the 
absolute error are given for the range. For example, 
absolute error is given for sin(jc) for values of x near tt. 

Sample: This column indicates the type of sample 
used for the accuracy figures. The type of sample de- 
pends upon the function and range under consid- 
eration. The statistics may be based either upon an 
exponentially (E) distributed argument sample or a 
uniformly (U) distributed argument sample.. 

Accuracy Figures: This column gives accuracy fig- 
ures for one or more representative segments within 
the valid argument range. The accuracy figures sup- 
plied are based upon the assumption that the argu- 
ments are perfect (i.e., without error and, therefore. 



having no error propagation effect upon the answers ) . 
The only error in the answers are those introduced by 
the subprograms. Appendix A contains a description 
of some of the symbols used in this appendix; the 
following additional symbols are used in the presenta- 
tion of accuracy figures: 



M (e) = Max 



i(x) - g(x) 



f(x) 



The maximum 
relative error 
produced 
during testing. 



(e) 



■4 






The standard 

f ( Xi ) - g(Xi) 2 deviation ( root- 
mean-square) of 
the relative error. 



f(Xi) 



M(E) = Max|f(x) -g(x) 



(E) 






N^i 



f(^i) - g(^i) 



The maximum 
absolute error 
produced 
during testing. 

The standard 
deviation (root- 
mean-square) of 
the absolute error. 



In the formulas for the standard deviation, IV repre- 
sents the total number of arguments in the sample; 
i is a subscript that varies from 1 to N. 

Accuracy for the Model 44 is based on performance 
with the FLOATING-POINT PRECISION switch in the "14" 
(vertical) position. This position selects the highest 
(56) of the four long-precision increments permitted. 

Average Speed: This column gives the timing statis- 
tics. These statistics represent the average speed in 
microseconds for the various System/ 360 models. Sta- 
tistics are supplied for Models. 30, 40, 44, 50, 65, and 
75. Statistics for the Model 75 are based upon two-way 
interleaving. 
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Table 12. Performance Statistics 
























1 

Entry 
Name 


Argument 
Range 


Sample 

E/U 


Accuracy Figures 


Average Speed ( Microseconds ) 


M(e) 


aie) 


M(E) 


o-(E) 


30 


40 


50 


65 1 75 
(See Note 8) 


44 
(Sec Note 9) 


ALGAMA 

ALOG 

ALOGIO 

ARCCTS 


< X ^ 0.5 


U 


1.09 X 10-6 


3.35 X 10-7 






10500 


2800 


865 


221 


131 


727 553 


0.5 < X <3 

3^x<« 
8^x< 16 
16^x<500 
0.5 ^ X ^ 1.5 
X < 0.5, X > 1.5 


u 

u ■■" 

~" u ' 

u 
" u" 






9.65 X 10-7 


3.74 X 10-7 


10500 


2820 


884 


225 


133 


734 560 


1.21 X 10-6 


2.86 X 10-7 






12100 


3250 


1020 


259 


151 


820 638 


1.25 X 10-6 


3.89 X 10-7 


3.46 X 10-7 




7600 


2010 


617 


162 


97.3 


470 356 


1.04 X 10-6 


2.03 X 10-7 




7600 


2010 


617 


162 


97.3 


477 362 






8.62 X 10-8 


4481 


1178 


361 


91.9 


52.3 


229 176 


E 


8.32 X 10-7 


1.20 X 10-7 


1.64X10-7^ 




4481 


1178 


361 


91.9 


52.3 


229 176 


0.5 ^ X ^ 1.5 
X < 0.5, X > 1.5 
-l^x:;^ + l 


U 
^ E " 

" U ' 

u~ 

Note 7 






4.78 X 10-8 


4847 


1278 


388 


98.1 


56.1 


248 193 


1.05 X 10-6 


2.17 X 10-7 




4847 


1278 


388 


98.1 


56.1 


248 193 


1.80 X 10-7 


3.29 X 10-6 


- - 




5340 


1460 


451 


118 


70.6 


325 238 


ARSIN 


- 1 g X =S + 1 


8.56 X 10-7 
9.75 X 10-7 


2.38 X 10-7 
4.54 X 10-7 

4.54 X 10-7 


5270 


1450 


445 


114 


68.4 


312 228 


ATAN 


The full range 
(Mod. 44: 
tan y, y in 

TT TT \ 
— '—-) 




For ATAN in FORTRAN IV (E) 
3602 913 255 68.8 40.6 
Additional time for ATAN in 

FORTJ^AN IV 
104 47 24 8 5.9 


165 125 


ATAN2 


The full range 


Note 7 


9.75 X 10-7 






4874 


1288 


375 


103 


19.7 




CABS 


The full range 


Note 1 


1.87 X 10-6 


7.65 X 10-7 


■ ■ - --■■ 




5194 


1421 


396 


108 


68.4 


299 205 


CCOS 


|x,| g 10, jx^l ^ 20 


U 


1.79 X 10-6 
See Note 2 


7.21 X 10-7 




15783 


4433 


1329 


334 


203 


1089 848 


CDABS 
CDCOS 


The full range 


Notel 


3.32 X 10-15 


5.16 X 10-16 







14021 


3061 


638 


150 


89.0 


733 656 


|xi| ^ 10, jx^l ^ 1 


U 


5.16 X 10-15 
See Note 3 


3.42 X 10-16 




48343 


11705 


2335 


507 


301 


3020 2807 


CDEXP 


Ix.| ^ 1, |x.| g 2" 


U 


4.04 X 10-16 


1.39 X 10-16 






41737 


10144 


2048 


456 


260 


2616 2419 


|xi| ^ 20, Ixal ^ 20 


U 


3.63 X 10-15 


1.29 X 10-16 






41737 


10144 


2048 


456 


260 


2611 2415 


CDLOG 


The full range 


Notel 


8.73 X 10-15 


6.38 X 10-17 


-— 




53362 


11940 


2393 


542 


317 


3105 2860 


CDSIN 


|xi| ^ 10, Ix^l ^ 1 


U 


3.72 X 10-15 
See Note 4 


3.49 X 10-16 




48250 


11668 


2322 


503 


298 


3018 2808 


CDSQRT 


The full range 


Notel 


9.86 X 10-16 


1.91 X 10-16 


- - - -- 


27951 


5906 


1282 


301 


181 


1543 1392 


CEXP 
CLOG 


|xi| gl70, |x2| ^j 


U 


1.18 X 10-6 


2.34 X 10-7 


13731 


3888 


1166 


291 


177 


958 732 


|xi| ^ 170, 

f <N:g20 


U 


1.06 X 10-6 


2.51 X 10-7 




13899 


3930 


1180 


294 


178 


971 746 


The full range 
except (1 + 01) 


Notel 


2.00 X 10-6 


1.56 X 10-7 






15504 


4246 


1261 


338 


216 


1014 777 



t^ Table 12. Performance Statistics ( Continued ) 



Entry 
Name 


Argument 
Range 


Sample 

E/U 


Accuracy Figures 


Average Speed ( Microseconds ) 


M(e) 


cr(e) 


M (E) 


a(E) 


30 


40 


50 


65 1 75 
(See Note 8) 


44 
(See Note 9) 


COS 


O^X^TT 


U 






1.47 X 10-7 


5.48 X 10-8 


3934 


1047 


298 


74.2 


44.0 


196 


157 


-10gx<0 

TT < X ^ 10 


u 






1.42 X 10-7 


5.67 X 10-8 


3990 


1061 


303 


75.4 


44.5 


198 


159 


10 < lx| ^ 100 


u 






1.35 X 10-7 


5.61 X 10-8 


3990 


1061 


303 


75.4 


44.5 


200 


161 


COSH 


-5^x^ +5 


u 


1.31 X 10-6 


3.40 X 10-6 






6110 


1810 


570 


145 


89.2 


486 


312 


COTAN 


1 1 <r- ■^ 

N^4 


u 


1.29 X 10-6 


3.68 X 10-7 






4420 


1180 


341 


86.7 


56.0 


227 


180 




JL<\x\^JL 
4 2 


u 


3.80 X 10-4 
See Note 5 


7.70 X 10-4 






4610 


1220 


351 


89.3 


55.5 


233 


188 


~<\x\ <10 
2 ' ' ~ 


u 


1.13 X 10-5 
See Note 5 


6.03 X 10-7 






4580 


1210 


348 


88.1 


55.0 


233 


188 


10 < |x| ^ 100 


u 


1.67 X 10-5 
See Note 5 


6.67 X 10-7 






4580 


1210 


348 


88.1 


55.0 


233 


187 


CSIN 


lxi| ^ 10, |X2| ^ 1 


u 


1.97 X 10-6 
See Note 6 


7.09 X 10-7 






15690 


4397 


1316 


331 


200 


1081 


843 


CSQRT 


The full range 


Notel 


1.61 X 10-6 


4.58 X 10-7 






10408 


2870 


805 


219 


140 


676 


493 


DARCOS 


-1^x^+1 


U 


2.72 X 10-16 


9.35 X 10-17 






22600 


5100 


1100 


246 


143 


1439 


1289 


DARSIN 


-1^x^+1 


U 


2.40 X 10-16 


6.00 X 10-17 






22400 


5060 


1090 


243 


140 


1440 


1292 


DATAN 


The full range 
(Mod. 44: 
tan y, y in 

TT TT \ 


Note 7 


2.08 X 10-16 


6.64 X 10-17 






For DATAN in FORTRAN IV (E) 

19056 4000 715 153 83.6 

Additional time for DATAN in 

FORTRAN IV 

104 47 24 8 5.9 


1010 


966 


DATAN2 


The full range 


Note 7 


2.08 X 10-16 


6.64 X 10-17 






22281 


4734 


886 


195 


22.9 








DCOS 


O^X^TT 


U 






1.79 X 10-16 


6.40 X 10-17 


13133 


3146 


605 


132 


74.2 


761 


713 


- 10 < X < 
7r<x^l0 


U 






1.76 X 10-16 


5.93 X 10-17 


13133 


3146 


605 


132 


74.2 


762 


714 


10 < |x| ^ 100 


U 






2.65 X 10-15 


1.01 X 10-15 


13133 


3146 


605 


132 


74.2 


759 


712 


DCOSH 


-5^x^ + 5 


u 


4.81 X 10-16 


1.34 X 10-16 






18300 


4260 


898 


206 


119 


1050 


923 


DCOTAN 


1 1 <- 'r 


u 


3.46 X 10-16 


8.38 X 10-17 






15300 


3590 


675 


150 


89.2 


869 


819 


-<|xl<- 
4 ^1^^1-2 


u 


1.72 X 10-13 
See Note 5 


5.00 X 10-15 






15900 


3640 


715 


156 


89.8 


906 


856 


^< |x| ^ 10 


u 


5.33 X 10-13 
See Note 5 


1.09 X 10-14 






15800 


3690 


706 


155 


89.1 


901 


852 


10 < |x| ^ 100 


u 


8.61 X 10-13 
See Note 5 


4.61 X 10-14 






15800 


3690 


706 


155 


89.1 


899 


849 



Table 12. Performance Statistics ( Continued ) 



> 



Entry 
Name 


Argument 
Range 


Samjile 

E/U 


Accuracy Figures 


Average Speed (Microseconds) 


M(e) 


(/(e) : 


M(E) 


(r(E) 


30 


40 


50 


65T^5 
(See Note 8) 


44 
( See Note 9 ) 


DERF 


|x| ^1.317 


U 


1.70 X 10- 1« 


2.71 X 10-17 






18400 


4610 


879 


190 


109 


1208 


1083 


1.317 < |x| ^ 2.04 


U 


2.91 X 10-^" 


1.17 X 10-17 






24300 


6060 


1150 


250 


141 


1635 


1467 


2.04 < |x < 6.092 


u 


1.70 X 10-17 


8.03 X 10-i« 






45200 


10900 


2080 


449 


259 


2868 


2626 


DERFC 


- 6 < X < 


u 


1.88 X 10-i« 


6.84 X 10-17 






36700 


8920 


1700 


369 


218 


2355 


2147 


0:;^x^ 1.317 


u 


3.52 X 10-i« 


7.62 X 10-17 






18600 


4650 


891 


193 


111 


1213 


1086 


1.317 < X ^ 2.04 


u 


4.45 X 10-i« 


1.27 X 10-16 






24100 


6000 


1 130 


244 


139 


1623 


1454 


2.04 < X < 4 


u 


4.02 X 10-i^> 


1.24 X 10-15 






45000 


10800 


2060 


444 


258 "^ 


2854 


2612 


4 -^ X < 13.3 


u 


5.02 X 10-3 5 


1.40 X 10-15 






45200 


10900 


2090 


451 


263 


2871 


2628 


DEXP 


|x| ^" 1 


u 


2.27 X 10-10 


7.49 X 10-17 






12145 


2907 


607 


138 


75.5 


720 


648 


1< |x| ^ 20 


u 


2.31 X 10-15 


8.69 X 10-16 






12145 


2907 


607 


138 


75.5 


716 


644 


20<|x| r:: 170 


u 


2.33 X 10-15 


9.33 X 10-16 






12145 


2907 


607 


138 


75.5 


715 


643 


DGAMMA 


< X < ] 


u 


2.18 X 10-ie 


7.93 X 10-17 


__.,. ,^._ _„ 




30700 


7500 


1400 


304 


175 


2058 


1853 


1 ^ X ^ 2 


u 


3.12 X 10-17 


8.45 X 10-i« 




28200 


7050 


1340 


292 


169 


1921 


1714 


2 < X ^ 4 


u 


1".69 X 10-15 


4.60 X 10-16 






30100 


7520 


1430 


312 


180 


2034 


1827 


4 < X < 8 


u 


2.85 X 10-15 


9.46 X 10-16 






33900 


8470 


1610 


352 


205 


2268 


2054 


8 ^ X < 16 


u 


6.42 X 10-15 


2.01 X 10-15 






40600 


9560 


1940 


434 


241 


2498 


2294 


16 ^ X < 57 


u 


6.2 X 10-14 


2.96 X 10-14 






40600 


9560 


1940 


434 


241 


2503 


2298 


DLGAMA 


< X ^ 0.5 


u 


4.11 X 10-16 


1.60 X 10-16 






46900 


11300 


2160 


471 


267 


3056 


2783 


0.5 < X < 3 


u 






2.86 X 10-16 


1.16 X 10-16 


44900 


11100 


2130 


466 


264 


2983 


2709 


3 :i X < 8 


u 


2.38 X 10-15 


3.99 X 10-16 






45900 


12200 


2330 


512 


292 


3228 


2947 


8 =i X < 16 


u 


3.36 X 10-16 


1.18 X 10-16 






28200 


6580 


1310 


296 


161 


1712 


1572 


16 gx<500 


u 


1.62 X 10-15 


2.43 X 10-16 






28200 


6580 


1310 


296 


161 


1725 


1585 


DLOG 


0.5 ^ X ^ 1.5 


u 






1.85 X 10-16 


7.29 X 10-17 


16044 


3769 


734 


161 


86.5 


929 


855 


x < 0.5, X > 1.5 


E 


3.31 X 10- 1« 


5.46 X 10-17 


8.23 X 10-17 




16041 


3765 


733 


161 


86.5 


929 


855 


DLOGIO 


0.5 =^ X ^; 1.5 


u 






3.09 X 10-17 


17149 


4048 


778 


171 


92.3 


997 


920 


x<0.5, X > 1.5 


E 


6.14 X 10-16 


9.96 X 10-17 






17147 


4044 


777 


170 


92.0 


996 


920 


DSIN 


1 1 «r "" 


u 


4.08 X 10-16 


4.85 X 10-17 


9.10 X 10-17 


2.17 X 10-17 


13145 


3148 


609 


133 


75.5 


762 


714 


f<|x|^10 


u 






1.64 X 10-16 


6.35 X 10-17 


13145 


3148 


609 


133 


75.5 


763 


715 


10 < |x|g 100 


u 






2.69 X 10-15 


1.03 X 10-15 


13145 


3148 


609 


133 


75.5 


761 


714 


DSINH 


|x| ^ 0.34657 


u 


2.10 X 10-16 


5.29 X 10-17 






8650 


2170 


408 


88.5 


52.4 


505 


457 


0.34657 < !x| ^ 5 


u 


3.59 X 10-16 


8.73 X 10-17 






18400 


4280 


901 


207 


119 


1049 


925 



[fe Table 12. Performance Statistics ( Continued ) 



Entry 
Name 


Argument 
Range 


Sample 

E/U 


Accuracy Figures 


Average Speed ( Microseconds ) 


M(£) 


aU) 


M (E) 


cr(E) 


30 

■ 


40 


50 


65 1 75 
(See Note 8) 


44 
(See Note 9) 


DSQRT 


The full range 


E 


1.08 X 10 -i« 


2.17 X 10-17 






8173 


1684 


355 


85.3 


49.2 


370 


334 


DTAX 


waf 


U 


5.25 X 10-16 


9.26 X 10-17 






15100 


3500 


647 


142 


84.1 


852 


806 
854 


■-<!x'< — 
4 ^ -^ - 2 


u 


1.67 X 10-12 
See Note 5 


3.69 X 10-14 






15700 


3660 


696 


151 


86.8 


902 


~ < |xt ^ 10 


u 


1.57 X 10-i:< 
See Note 5 


4.51 X 10-13 






15600 


3640 


688 


150 


86.1 


895 


848 


10 < |x] ^ 100 


u 


3.79 X 10-1-' 
Sec Note 5 


9.50 X 10-14 






15600 


3640 


688 


150 


86.1 


893 


846 


DTANH 


|x| ^ 0.54931 


u 


2.00 X 10-1" 


4.45 X 10-17 






12299 


2850 


477 


106 


55.5 


668 


641 


0.54931 < |x| ^ 5 


u 


1.99 X 10-'"' 


2.54 X 10-17 






16078 


3778 


833 


192 


110 


979 


874 


EXP 


|x| ^ 1 


u 


4.65 X 10-7 


1.28 X 10-7 
1.17 X 10-7 






4173 


1250 


388 


95.2 


53.4 


311 


192 


|x| ^ 170 


u 


4.69 X 10-7 






4183 


1250 


387 


94.8 


53.4 


308 


189 


ERF 


lx| g 1.317 


u 


9.26 X 10-7 


1.43 X 10-7 






4140 


1120 


363 


88.5 


52.7 


276 


189 


1.317 <|x| ^2.04 


u 


9.02 X 10-« 


3.42 X 10 -« 






4500 


1230 


397 


99.1 


58.8 


322 


220 


2.04 < |x| ^ 3.9192 


u 
u 


6.07 X 10 -« 
9.10 X 10-7 


3.42 X 10-^ 
2.97 X 10-7 






10000 


2800 


845 


213 


127 


740 


525 


ERFC 


- 3.8 < x < 






7040 


1960 


607 


153 


91.6 


521 


367 


O^x;^ 1.317 


u 


3.90 X 10-" 


5.65 X 10-7 






4250 


1150 


374 


92.0 


54.3 


284 


195 


1.317 <x^ 2.04 ^ 


u 


1.02 X 10-" 
1.20 X 10-" 


2.13 X JO-7 






4410 


1210 


387 


96.0 


58.0 


319 


216 


2.04 < X < 4 


u 


3.60 X 10-7 






9950 


2780 


835 


210 


126 


737 


521 
"532 


4 ^ x ^ 13.3 


u 


1.52 X 10-"> 


8.45 X 10-" 






10100 


2840 


859 


216 


131 


749 


GAMMA 


< X < 1 


u 


9.86 X 10-7 


3.45 X 10-7 






5840 


1560 


484 


123 


74.0 


433 


306 


1 ^ X ^ 2 


u 


1.00 X 10-7 


3.74 X 10- s 






5620 


1520 


489 


123 


73.5 


428 


301 
330 


2 < X ^ 4 


u 


9.29 X JO-7 


3.63 X 10-7 






6330 


1700 


546 


137 


81.1 


460 


4 < X < 8 


u 


2.25 X 10-" 


8.14 X 10-7 






7740 


2070 


659 


166 


96.2 


538 


402 


8 ^ X < 16 


u 


2.29 X 10-"' 


7.67 X 10-" 






12000 


3320 


1020 


263 


155 


845 


618 


16 g X < 57 


u 


4.36 X 10—-' 


1.45 X 10- •' 






12000 


3320 


1020 


263 


155 


842 


616 


SIN 


1 1 «- "■ 


u 


1.59 X 10-" 


2.02 X 10-7 


1.31 X 10-7 


5.55 X 10-8 


3876 


1036 


298 


74.2 


44.3 


194 


155 


|-<|x|^10 


u 






1.41 X 10-7 


5.53 X 10-« 


3989 


1064 


307 


76.4 


45.3 


201 


163 


10 < |x| g 100 


u 






1.46 X 10-7 


5.61 X 10- s 


3989 


1064 


307 


76.4 


45.3 


200 


162 


SINH 


- 5^x^ + 5 


u 


1.20 X 10-" 


3.20 X 10-7 






5890 


1740 


545 


139 


85.3 


461 


297 


SQRT 


The full range 


E 


8.70 X 10-7 


1.68 X 10-7 






2965 


801 


210 


59.1 


35.8 


142 


90 



Table 12. Performance Statistics ( Continued ) 



> 



Entry 
Name 


Argument 
Range 


Sample 

E/U 




Accuracy Figures 




Average Speed ( Microseconds ) 


M(c) 


<rie) 


M (E) 


C7(E) 


30 


40 


50 


65 1 75 
(See Note 8) 


44 
( Sec Note 9 ) 


TAN 
TANK 


wsj 


U 


1.56 X 10-6 


3.22 X lO-T 


— 




4220 


1120 


319 


79.9 


51.3 


209 166 


4 ^ '"' = 2 


U 


6.58 X 10-5 
See Note 5 


1.67 X 10-6 




4500 


1184 


338 


85.3 


52.8 


232 188 


f <Waio 


u 


4.92 X 10-5 
See Note 5 


1.28 X 10-6 




4460 


1170 


335 


84.1 


52.4 


226 183 


10 < |x| g 100 
|x| ^ 0.54931 


u 


3.35 X 10-5 
See Note 5 


1.02 X 10-6 




4460 


1170 


335 


84.1 


52.4 


227 184 


u 


8.12 X 10-7 


1.66 X 10-1 




2581 


649 


173 


46.1 


29.0 


89 63 


0.54931 < |x| ^ 5 


u 


5.74 X 10-7 


7.53 X 10-8 




5952 


1774 


551 


142 


86.3 


446 294 



Notes to Table 12 

These notes are associated with Table 12 and contain 
more detailed information about samples and relative 
errors for certain functions. 

Note 1: The distribution of sample arguments upon 
which these statistics are based is exponential radially 
and is uniform around the origin. 

Note 2: The maximum relative error cited for the 
ccos function is based upon a set of 2000 random argu- 
ments within the range. In the immediate proximity of 

1_ 

2 

. . . ,) the relative error can be quite high, although 
the absolute error is small. 



the points 



(-i)' 



+ Oi (where n= 0, ±1, ±2, 



Note 3: The maximum relative error cited for the 
CDCOs function is based upon a set of 1500 random 
arguments within the range. In the immediate prox- 
imity of the points I n -H -^ j tt + Ot ( where n == 0, 

±1, ± 2, . . . , ) the relative error can be quite high 
although the absolute error is small. 



Note 4: The maximum relative error cited for the 
CDsiN function is based upon a set of 1500 random 
arguments within the range. In the imnaediate prox- 
imity of the points nw + Oi (where n = ±1, ±2, 
. . . , ) the relative error can be quite high although 
the absolute error is small. 

Note 5: The figures cited as the maximum relative 
errors are those encountered in a sample of 2500 ran- 
dom arguments within the respective ranges. See the 
appropriate section in Appendix A for a description 
of the behavior of errors when the argument is near a 
singularity or a zero of the function. 

Note 6: The maximum relative error cited for the 
csiN function is based upon a set of 2000 randon argu- 
ments within the range. In the immediate proximity 
of the points n-n- + 0/ (where n = ±1, ±2, . . . ,) the 
relative error can be quite high although the absolute 
error is small. 

Note 7: The sample arguments were tangents of 
numbers uniformly distributed between — — and + -^ . 

Note 8: The statistics for the Model 75 are based 
upon two-way interleaving. 

Note 9: The second column of speeds for the Model 
44 applies to that machine with high-speed registers. 



Appendix C. Interruption and Error Procedures 



Appendix C contains descriptions of the procedures 
followed when the execution of a load module is dis- 
continued. Execution may be discontinued due to one 
of two reasons: an interruption or an error. After an 
interruption is processed, execution of this load module 
or phase continues; after an error is processed, execu- 
tion of this load module or phase is terminated. The 
following text explains the procedure used to handle 
each case. 



Interruption Procedures 

An interruption is a computer-originated break in the 
flow of processing. When an interrupt occurs, an indi- 
cator is set to record exponent overflow, underflow, or 
divide exception. ( The status of the indicators can be 
tested by using the subprograms described in "Service 






IKJ/i^LCLXXM 



in«.crrup«, message is men 



written. After the interruption is handled, execution of 
the load module ( or phase, in the Model 44 Program- 
ming System) continues from the point at which it 
was interrupted. 

The program interrupt message contains the old pro- 
gram status word (psw), which indicates the cause 
of the interrupt. Figure 1 shows the format of the 
message as it is issued by the operating system. Figure 
2 shows the format as issued by the Model 44 system. If 
the letter A appears in parentheses in the program 
interrupt message, boundary adjustment has taken 
place. The eighth character in the psw represents the 
code number associated with the type of interruption. 
These interruptions are described in the following 
paragraphs. (For more information on the psw, see 
the publication IBM System/ 360 Principles of Opera- 
tion, Form A22-6621. ) 

Specification Exception (Code 6): The specification 
exception ( code 6 ) is recognized when a data address 



does not specify an integral boundary for that unit of 
information. A specification error would occur, for ex- 
ample, during the execution of the following program 
segment: 

DOUBLE-PRECISION D, E 
COMMON A, B, C 
EQUIVALENCE (B, D) 
D = 3.0D02 

Fixed-Point Divide Exception (Code 9): The fixed- 
point divide exception (code 9) is recognized when 
division of a fixed-point number by zero is attempted. 
A fixed-point divide exception would occur during exe- 
cution of the following statements : 

J =0 
I = 7 
K = I/J 

Exponent-Overflow Exception (Code C): The ex- 
ponent-overflow exception, (code C) is recognized 
when the absolute value of the result of a floating- 
point addition, subtraction, multiplication, or division 
is greater than or equal to 16*^^ (approximately 7.2 x 
10^^ ) . For example, an exponent overflow would occur 
during execution of the statement: 

A = l.OE -I- 75 + 7.2E + 75 

Exponent-Underflow Exception (Code D): The ex- 
ponent-underflow exception, (code D) is recognized 
when the absolute value of the result of a floating-point 
addition, subtraction, multiplication, or division, is less 
than lO"**"' (approximately 5.4 x 10~"^) but not equal 
to 0. An exponent-underflow exception would occur 
during execution of the statement: 



A = - LOE - 50 



l.OE - 50 



Floating-Point Divide Exception (Code F): The 
floating-point divide exception (code F) is recognized 
when division of a floating-point number by zero is 



IHC210I PROGRAM INTERRUPT (A) OLD PSW IS xxxxxxx 



• Figure 1. Format of Program Interrupt Message, Operating System 
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OA210I PROGRAM INTERRUPT (A) OLD PSW IS xxxxxxx 



B 



xxxxxxxx 



Figure 2. Format of Program Interrupt Message, Model 44 System 

attempted. A floating-point divide exception would 
occur during execution of the following statements: 

B = 0.0 
A = 1.0 
C = A/B 



System /360 Operating System 

(The following paragraphs do not apply to the 
Model 44 Programming System. ) 

The interrupt message is written in the system out- 
put data set. 

A specification exception program interrupt message 
(code 6) is issued only if the boundry = align option 
was specified in the fortlib macro-instruction during 
system generation and a boundary afignment error 
occurs. Then the boundary alignment routine is in- 
voked to correct the boundary misalignment. If an 
instruction that has been processed for boundary mis- 
alignment also contains a protection, addressing, or 
data error, the interrupt message will be reissued with 
the appropriate code ( 4, 5, or 7 ) . (In these cases, the 
letter A appears in parentheses in the program inter- 
rupt message.) Then the job will terminate because 
both a specification error and a protection, addressing, 
or data error have been detected. The completion code 
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cause of the specification error. 

The number of warning messages printed is limited 
to ten. After ten boundary alignment adjustments have 
been made, the message is suppressed, but boundary 
alignment violations continue to be corrected. 

Protection Exception (Code 4): The protection ex- 
ception (code 4) is recognized when the key of an 
operand in storage does not match the protection key 
in the psw. A message is issued only if a specification 
exception ( code 6) has already been recognized in the 
same instruction. Otherwise, the job terminates abnor- 
mally without a message. 

Addressing Exception (Code 5): The addressing ex- 
ception (code 5) is recognized when the address of 
the data is outside of available storage for the particu- 
lar installation. A message is issued only if a specifica- 
tion exception (code 6) has already been recognized 
in the same instruction. Otherwise, the job terminates 
abnormally without a message. 



Data Exception (Code 7): The data exception ( code 
7) is recognized v^'hen the sign or digit codes for a 
CONVERT TO BINARY instruction are incorrect. A message 
is issued only if a specification exception (code 6) 
has already been recognized in the same instruction. 
Otherwise, the job terminates abnormally without a 
message. 

Model 44 Programming System 

( The following paragraphs do not apply to System/ 
360 Operating System. ) 

The interrupt message is written in sysopt. 

The program interrupt message (with code 6, as 
described in "Specification Exception") that results 
when the boundary specification convention is vio- 
lated will contain the "(A)" only if the &fix option 
has been turned on (seta i) in boauopt. With that 
option on ( as it is in the distributed version of the sys- 
tem), a routine that adjusts for the misalignment is 
executed each time such a violation occurs, and proc- 
essing continues. 

With the &FIX option on, the number of program 
interrupt messages put out as the result of boundary 
violations is limited to a message for each of the first 
n violations per execution, where n is the operand of 
the SETA instruction for the &prntmes option. In the 
system as distributed, this operand is equal to 0. Only 
the message, and not the alignment correction, is 
inhibited. 

Errot Procedures 

During execution, the mathematical subprograms as- 
sume that the argument(s) is the correct type. No 
checking is done for erroneous arguments (i.e., the 
wrong type, invalid characters, the wrong length, etc. ) ; 
therefore, a computation performed with an erroneous 
argument has an unpredictable result. However, the 
nature of some mathematical functions requires that 
the input be within a certain range. For example, the 
square root of a negative number is not permitted. 
If the argument is not within the valid range given in 
Tables 2 through 6, an error message is written in the 
data set associated with system output. The execution 
of this load module or phase is terminated and control 
is returned to the operating system. 
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The error message that is issued has the format : 

IHCyyyl 

TRACEBACK FOLLOWS . . . 

or 
OAyyyl 

The former message is issued by the operating system, 
the latter by the Model 44 system. Traceback is a diag- 
nostic tool for the Operating System Fortran iv 
only. It is a list of routines in the direct line of call to 
the routine in which the error occurred. It is described 
in either IBM System/360 Operating System FOR- 
TRAN IV Programmers Guide: FORTRAN IV (G), 
Form No. C28-6639, or FORTRAN IV (H), Form No. 
C28-6692. The yyy is a numeric code that identifies the 
error detected. The following text lists the error mes- 
sages in numeric order, explains the error, and indi- 
cates what action the system takes. In the following 
explanations, x represents the argument supplied by 
the programmer. 

IHC216l;OA216l 

Explantation: In the xxxfslit subprogram, a value of 
i that is not 0, 1, 2, 3, or 4 is an error. 
System Action: Execution of this load module or 
phase is terminated. 

IHC241l;OA241l 

Explanation: In the xxxfixpi subprogram, a base 
number of zero and an exponent ^ is an error. 
System Action: Execution of this load module or 
phase is terminated. 

IHC242I; OA242I 

Explanation: In the xxxfrxpi subprogram, a base 
number of zero and an exponent ^ is an error. 
System Action: Execution of this load module or 
phase is terminated. 

IHC243I; OA243I 

Explanation: In the ocxxfdxpi subprogram, a base 
number of zero and an exponent ^ is an error. 
System Action: Execution of this load module or 
phase is terminated. 

IHC244I; OA244I 

Explanation: In the xxxfrxpr subprogram, a base 
number of zero and an exponent ^ is an error. 
System Action: Execution of this load module or 
phase is terminated. 

IHC2451; OA245I 

Explanation: In the arxxFoxPD subprogram, a base 
number of zero and an exponent ^ is an error. 
System Action: Execution of this load module or 
phase is terminated. 



IHC246I; OA246I 

Explanation: In the xxxfcxpi subprogram, a base 
number of zero and an exponent ^ is an error. 
System Action: Execution of this load module or 
phase is terminated. 

IHC247I; OA247I 

Explanation: In the xxxfcdxi subprogram, a base 
number of zero and an exponent ^ is an error. 
System Action: Execution of this load module or 
phase is terminated. 

IHC251l;OA251l 

Explanation: In the xxxssqrt subprogram, a value 
of X < is an error. 

System Action: Execution of this load module or 
phase is terminated. 

IHC252I; OA252I 

Explanation: In the xxxsexp subprogram, a value of 
X > 174.673 is an error. 

System Action: Execution of this load module or 
phase is terminated. 

IHC2S3I; OA253I 

Explanation: In the jcxxslog subprogram, a value of 
X ^ is an error. Because this subprogram is also 
called by an exponentiation subprogram, this mes- 
sage also indicates that an attempt has been made 
to raise a negative real base to a power. 
System Action: Execution of this load module or 
phase is terminated, 

IHC25,4h OA254I 

Explanation: In the xxxsscn subprogram, a value of 
!x!^ 2^^ • TT is an error. 

System Action: Execution of this load module or 
phase is terminated. 

IHC255I; OA255I 

Explanation: In the xxxsatn2 subprogram when 
entry name atan2 is used, a value of x, = Xo = is 
an error. 

System Action: Execution of this load module or 
phase is terminated. 

IHC256h 0'A2S6I 

Explanation: In the xxxsscnh subprogram, a value 
of [x| ^ 174.673 is an error. 

System Action: Execution of this load module or 
phase is terminated. 

IHC257I! OA257I 

Explanation: In the xxxsascn subprogram, a value 
of |x[ > 1 is an error. 

System Action: Execution of this load module or 
phase is terminated. 
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IHC258I; OA258I 

Explanation: In the xxxstnct subprogram, a value of 
Lt| ^ 2''"' • 77 is an error. 

System Action: Execution of this load module or 
phase is terminated. 

IHC259I; OA259I 

Explanation: In the xxxstnct subprogram, a value of 

X too close to one of the singularities ( ± -^, ± -r-, . . . 

for the tangent; ±7r, ±27r, . . . for the cotangent) is 
an error. 

System Action: Execution of this load module or 
phase is terminated. 

IHC26U;OA261l 

Explanation: In the xxxlsqrt subprogram, a value 
of X < is an error. 

System Action: Execution of this load module or 
phase is terminated. 

IHC262I; OA262I 

Explanation: In the xxxlexp subprogram, a value of 
X > 174.673 is an error. 

System Action: Execution of this load module or 
phase is terminated. 

IHC263I; OA263I 

Explanation: In the xxxllog subprogram, a value of 
X ^ is an error. Because this subprogram is also 
called by an exponentiation subprogram, this mes- 
sage also indicates that an attempt has been made 
to raise a negative real number to a power. 
System Action: Execution of this load module or 
phase is terminated. 

IHC264I; OA264I 

Explanation: In the xxxlscn subprogram, a value of 
'xj ^ 2'*" • 7- is an error. 

System Action: Execution of this load module or 
phase is terminated. 

IHC265I; OA265I 

Explanation: In the xxxlatns subprogram when 
entry name datan2 is used, a value of Xi = X2 = 
is an error. 

System Action: Execution of this load module or 
phase is terminated. 

IHC266I; OA266I 

Explanation: In the xxxlscxh subprogram, a value 
of |x; ^ 174.673 is an error. 

System Action: Execution of this load module or 
phase is terminated. 



IHC2671; OA267I 

Explanation: In the xxxlascn subprogram, a value 
of jx! > 1 is an error. 

System Action: Execution of this load module or 
phase is terminated. 

IHC268I! OA268I 

Explanation: In the xxxltnct subprogram, a value of 
ix! ^ 2^° • 77 is an error. 

System Action: Execution of this load module or 
phase is terminated. 

IHC269I; OA269I 

Explanation: In the xxxltnct subprogram, a value of 

X too close to one of the singularities {—~k', —~^> • • • 

for the tangent; ±77, ± £77, . . . for the cotangent ) is 
an error. 

System Action: Execution of this load module or 
phase is terminated. 

IHC271l;OA271l 

Explanation: In the xxxcsexp subprogram, a value 
of Xi > 174.673 is an error. 

System Action: Execution of this load module or 
phase is terminated. 

mC272l; OA272I 

Explanation: In the xxxcsexp subprogram, a value 
of |x2| ^ 2^^ • 77 is an error. 

System Action: Execution of this load module or 
phase is terminated. 

IHC273I; OA273I 

Explanation: In the xxxcslog subprogram, a value 
Oj. X]i -— X2 "~" v/ IS an error. 

System Action: Execution of this load module or 
phase is terminated. 

IHC2741; OA274I 

Explanation: In the xxxcsscx subprogram, a value 
of [xij ^ 2^^ • 77 is an error. 

System Action: Execution of this load module or 
phase is terminated. 

IHC275I; OA275I 

Explanation: In the xxxcsscn subprogram, a value 
of |x2|> 174.673 is an error. 

System Action: Execution of this load module or 
phase is terminated. 

/HC28I/;OA28I/ 

Explanation: In the xxxclexp subprogram, a value 
of X] > 174.673 is an error. 

System Action: Execution of this load module or 
phase is terminated. 
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IHC282I! OA282I 

Explanation: In the xxa:cLEXP subprogram, a value 
of 1x2 1 ^ 2'^° • TT is an error. 

System Action: Execution of this load module or 
phase is terminated. 

IHC283I; OA283I 

Explanation: In the xxxcllog subprogram, a value 
of Xi = X2 = is an error. 

System Action: Execution of this load module or 
phase is terminated. 

IHC284I; OA284I 

Explanation: In the xxxclscn subprogram, a value 
of IxJ ^ 2^^ • TT is aA error. 

System, Action: Execution of this load module or 
phase is terminated. 

IHC285I; OA285I 

Explanation: In the xxxclscn subprogram, a value 
of |x2| > 174.673 is an error. 

phase is terminated. 



IHC2901; OA290I 

Explanation: In the xxxsgama subprogram for the 
gamma function, a value of x ^ 2--^- orx ^ 57.5744 
is an error. 

System Action: Execution of this load module or 
phase is terminated. 

/HC29I/,OA29?/ 

Explanation: In the xxxsgama subprogram for the 
log-gamma function, a value of x ^ or x ^ 4.2937 



10^ 



IS an error. 



System Action: Execution of this load module or 
phase is terminated. 

IHC300h OA300I 

Explanation: In the xxxlgama subprogram for the 
gamma function, a value of x ^ 2~^^- or x ^ 57.5744 
is an error. 

System Action: Execution of this load module or 
phase is terminated. 

IHC301I;OA301I 

Explanation: In the xxxlgama subprogram for the 
log-gamma function, a value of x ^ or x ^ 4.2937 
• 10'^'^ is an error. 

phase is terminated. 
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Appendix D. Storage Estimates 



Appendix D contains decimal storage estimates (in 
bytes ) for the library subprograms. The estimate given 
does not include any additional library subprograms 
or FORTRAN exccution-time routines that the subpro- 
gram needs during execution. The names of any addi- 
tional library subprograms needed are given in iabies 
13 and 14 in the column headed "Additional Subpro- 
grams." 

Some library subprograms also require execution- 
time routines for input/ output, interruption, and error 
procedures. 

If the programmer has not made allowances for the 
storage required by any of these additional routines 
( see Tables 15 and 16), the amount of available storage 
may be exceeded and execution cannot begin. The pro- 
grammer must add the estimates for all subprograms 
and routines needed to determine the amount of stor- 
age required. 



System/ 360 Operating System 

The iHCFiosH routine performs input/ output pro- 
cedures for both FORTRAN IV (e) and Fortran iv. [This 
routine refers to a table (ihcuatbl) for information 
about the input/ output devices used during execu- 
tion.] The IHCFCOME routine performs interruption and 
error procedures for Fortran rv (e) library subpro- 
grams; the IHCFCOMH, iHCFcvTH, iHCTRCH, and IHCUOPT' 
routines perform the procedures for Fortran rv library 
subprograms. If a system contains both compilers, the 
IHCFCOMH-IHCFCVTH routines are used. Tables 13 and 
14 indicate which library subprograms require these 
execution-time routines. 

In addition, several other execution-time routines 
may be needed to resolve external references in a 
FORTRAN IV object modulc. 

1. If a source module specifies direct-access input/ 
output operations, the compiler generates a call to 
the iHCDiosE routine. 

2. At the point that errors are encountered during 
compilation, the compiler generates a call to an 
error routine (ihciberr for Fortran iv (e) and 
iHciBERH for FORTRAN iv). If exccution of the load 
module is attempted, the error routine is called, 
a message is issued, and execution is terminated. 

3. If a FORTRAN IV (e) source module contains a com- 
puted GO TO, the compiler generates a call to the 
HCCGOTO routine. 



Table 13. Mathematical Subprogram Storage Estimates 











Uses 










Input/ 




OS/360 


44 PS 




Output 




Deci- 


Deci- 




and 




mal 


mal 




Inter- 


Subprogram 


Esti- 


Esti- 


Additional 


ruption 


Name 


mate 


mate 


Subprograms 


Routines 


xxxCLABS 


170 


200 


xxxLSQRT 


Yes 


xxxCLAS 


210 


220 




No 


xxxCLEXP 


250 


280 


xxxLEXP, xxxLSCN 


Yes 


xxxCLLOG 


260 


310 


xxxCLABS, xxxLSQRT, 
xxxLLOG, xxxLATN2 


Yes 


xxxCLSCN 


400 


500 


xxxLEXP, xxxLSCN 


Yes 


xxxCLSQT 


200 


240 


xxxLSQRT 


Yes 


xxxCSABS 


160 


190 


xxxSSQRT 


Yes 


xxxCSAS 


190 


200 




No 


xxxCSEXP 


240 


280 


xxxSEXP, xxxSSCN 


Yes 


xxxCSLOG 


240 


290 


xxxSABS, xxxSSQRT 
xxxSLOG, XXXSATN2 


Yes 


xxxCSSCN 


380 


440 


xxxSEXP, xxxSSCN 


Yes 


xxxCSSQT 


190 


230 


xxxSSQRT 


Yes 


IHCFAINT 


80 


— 




No 


xxxFCDXI 


300 


370 


xxxCLAS 


Yes 


xxxFCXPI 


280 


350 


xxxCSAS 


Yes 


xxxFDXPD 


210 


240 


xxxLLOG, xxxLEXP 


Yes 


xxxFDXPI 


160 


200 




Yes 


IHCFIFIX 


120 







No 


xxxFIXPI 


170 


230 




Yes 


xxxFMAXD 


110 


170 




No 


xxxFMAXI 


210 


290 




No 


xxxFMAXR 


210 


290 




No 


IHCFMODR 


120 







No 


IHCFMODI 


60 







No 


xxxFRXPI 


150 


190 




Yes 


xxxFRXPR 


210 


250 


xxxSLOG, xxxSEXP 


Yes 


xxxLASCN 


400 


520 


xxxLSQRT 


Yes 


IHCLATAN 


320 







No 


XXXLATN2 


500 


520 




Yes 


xxxLERF 


800 


920 


xxxLEXP 


Yes 


xxxLEXP 


460 


480 




Yes 


xxxLGAMA 


730 


820 


xxxLLOG, xxxLEXP 


Yes 


xxxLLOG 


380 


430 




Yes 


xxxLSCN 


380 


390 




Yes 


xxxLSCNH 


230 


400 


xxxLEXP 


Yes 


xxxLSQRT 


150 


160 




Yes 


xxxLTANH 


340 


350 


xxxLEXP 


Yes 


xxxLTNCT 


390 


400 




Yes 


xxxSASCN 


300 


380 


xxxSSQRT 


Yes 


IHCSATAN 


200 







No 


XXXSATN2 


360 


380 




Yes 


xxxSERF 


450 


560 


xxxSEXP 


Yes 


xxxSEXP 


290 


340 




Yes 


xxxSGAMA 


510 


600 


xxxSLOG, xxxSEXP 


Yes 


xxxSLOG 


270 


280 




Yes 


xxxSSCN 


260 


270 




Yes 


xxxSSCNH 


280 


340 


xxxSEXP 


Yes 


xxxSSQRT 


180 


180 




Yes 


xxxSTANH 


270 


280 


xxxSEXP 


Yes 


xxxSTNCT 


290 


310 




Yes 
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4. If a FORTRAN IV source module contains any input/ 
output operations that refer to a namelist name, 
compiler generates a call to the ihcnamel routine. 

5. If a FORTRAN IV sourcc module uses the debug 
facility, the compiler generates a call to the ihcdbug 
routine- 

6. If boundary alignment was specified during system 
generation, the ihcadjst routine v^^ill be loaded if a 
boundary-alignment error occurs. 

Model 44 Programming System 

In the FORTRAN library of the Model 44 Programming 
System, the boafiocs routine is the interface with the 
system input/ output services. This routine refers to a 
table, BOAUNrrB, for information about the input/ output 
devices used during ekecution. The boaibcom routine 
performs interruption and error procedures. If a source 
program contains any input/ output operation(s) re- 
ferring to a NAMELIST name, the compiler generates a 
call to the boanamel routine. 

Table 14. Service Subprogram Storage Estimates 



Table 15. Execution-Time Routine Storage Estimates, Oper- 
ating System 





OS 


44 PS 


Uses Input/Output 


Ci.V>*-.*>r.#*.-/>»v> 


r»^^;w,„1 


T^ ,-^^1 


^^J T«4-««»«..^4-i^« 










Name 


Estimate 


Estimate 


Routines 


xxxFDVCH 


80 


120 


Yes 


xxxFDUMP 


450 


760 


Yes 


xxxFEXIT 


30 


40 


Yes 


xxxFOVER 


90 


130 


Yes 


xxxFSLIT 


190 


280 


Yes 



Routine 


Decimal 




Name 


Estimate 


Used By 


IHCADJST 


1,090 


FORTRAN IV 


IHCCGOTO 


60 


FORTRAN IV (E) 


IHCDBUG 


2,600 


FORTRAN IV 


IHCDIOSE 


2,500 


Both 


IHCFCOME 


5,500 


FORTRAN IV (E) 


IHCFCOMH 


4,050 


FORTRAN IV 


IHCFCVTH 


4,090 


FORTRAN IV 


IHCFIOSH 


3,800 + IHCUATBL 

(See Note) 


Both 


IHCIBERH 


210 


FORTRAN IV 


IHCIBERR 


260 


FORTRAN IV (E) 


IHCNAMEL 


2,250 


FORTRAN IV 


IHCTRCH 


590 


FORTRAN IV 


IHCUOPT 


8 


FORTRAN IV 


Note: The nu 


mber of bytes in table IHCUATBL may be 


computed by the formula 






12n-F 8 




where n 


is the number of data set reference numbers 


requested during system generation. 



Table 16. Execution-Time Routine Storage Estimates, Model 
44 System 



Routine Name 


Decimal Estimate 


BOADIOCS 


680 


jo^jAr iuc;a 


i,^oo 


BOAIBCOM 


9,172 


BOANAMEL 


2,520 


BOAUNITB 


( see Note 1 ) 


BOAUOPT 


8 


BNXADJST 


1,000 




(Note 2) 


Note 1: Th6 number of bytes in CSECT BOAUNITB may 


be computed by the formula 


8n -1-8 


where n is the number of data set reference numbers 


requested during system construction. 


Note 2: BNXADJST does not reside in the library of relo- 


catable modules, but in the absolute phase library. 
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Appendix E. Assembler Language Information 



The mathematical and service subprograms in the 
FORTRAN IV Hbrarv are available to the assembler lan- 
guage programmer. The following text explains the 
method of calling a library subprogram in an assembler 
language program, and then gives additional informa- 
tion necessary to use each type of subprogram. (The 
assembler language programmer should also be fami- 
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Calling Sequences 

To call either type of library subprogram, the assem- 
bler language programmer supplies an entry name, 
an argument list, and an area used by the subprogram 
to store information (i.e., a save area). The following 
conventions must be observed when calling a library 
subprogram in an assembler language program: 

1. The address of the entry name must be in general 
register 15. 

2. The address of the point of return to the calling 
program must be in general register 14. 

3. The address of the argument hst must be in general 
register 1. 

4. The argument list must be assembled on a full-word 
boundary; it consists of one 4-byte address constant 
for each argument. The last argument must have 
a 1 in its high order bit. 

5. The address of the save area must be in general 
register 13. 

6. The save area must be assembled on a full-word 
boundary. Although the minimum size of the save 
area depends upon the subprogram, the program- 
mer is advised to use a save area of 18 full-words 
for all library subprograms. The minimum save 
area sizes are given in Tables 2 through 6 for the 
mathematical subprograms, and in Table 17 for the 
service subprograms. 

7. If the information in a floating-point register is to 
be retained, the programmer must save and restore 
the contents of the register. The subprograms that 
make use of the floating-point registers contain no 
provisions for saving the information. 

8. If a main program in assembler language contains 
any calls to those library subprograms that use the 
FORTRAN execution-time routines ( see Appendix D ) , 
the following instructions must be included before 
the call to the subprogram is issued: 

OS 44PS 



L 15,=V(IBC0M#) 

BAL 14,64(15) 



EXTRX IBCOM# 



These instructions cause the initialization of return 
coding and the interruption exceptions described in 
Appendix C. If these instructions are omitted, the 
occurrence of an interruption or an error causes un- 
predictable termination of the execution of this load 
module. 
Note : In an assembler language program, a decimal- 
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acter B to appear in the program interruption 
message described in Appendix C. 

The user of System/ 360 Operating System may use 
several methods to call a Fortran library subprogram: 
the appropriate macro-instructions described in the 
publication IBM System/ 360 Operating System: Su- 
pervisor and Data Management Macro-Instructions, 
Form C28-6647 or the general assembler language 
calhng sequence (given in Figure 3). If the macro- 
instructions are used, the address of the save area must 
be placed in general register 13 before using a macro- 
instruction to give control to the subprogram. For 
example, if the square root of the value in amnt is to 
be computed and save is the address of the same area, 
the following statements could be included in an as- 
sembler language program to call the ihcssqrt sub- 
program: 

L 15,=V(IBC0M#) 

BAL 14,64(15) 



LA 
CALL 



13,SAVE 
SQRT,(AMNT),VL 



SAVE DS 18F 

If the general assembler language calling sequence 
shown in Figure 3 is used, the programmer must en- 
sure that all of the conventions discussed previously 
are followed. For example, to call the ihcssqrt sub- 
program to compute the square root of the number in 
amnt, the following statements would be included in 
the source program: 

LA 13,SAVE 

LA 1,ARG 

L 15,ENTRY 

BALR 14,15 



L 15, = A(IBCOM#) 

BAL 14,64(15) 



ENTRY 



SAVE 



ARC 



DC 



DS 



DC 
DC 



V(SQRT) 



18F 



X'80' 
AL3(AMNT) 
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entry 
entry 



arglist 



LA 13, area 

LA 1, arglist 

L 15, entry 

BALR 14, 15 



NOP 



X'id' 



DC V { entry name ) 

or 
DC A ( entry name ) 



DS 



xxF 



CNOP 

DC 



X'80' 



DC AL3 (arg) 

or for more than one argument: 
arglist DC A ( argi ) 



DC 



A (arga) 



General register 13 contains the address of the save area. 

General register 1 contains the address of the argument list. 

General register 15 contains the address of the subprogram. 

General register 14 contains the address of the point of return to the 
calling program. 

This statement is optional. The id represents the binary calling sequence 
identifier. This number is supplied by the programmer and may be any 
hexadecimal integer less than £16 — 1. 



Note : In this case, the entry name must be defined by an EXTRN instruc- 
tion to obtain proper linkage. 

This statement defines the save area needed by the subprogram. The xx 
represents the minimum size of the save area required; however, the pro- 
grammer is advised to use a save area of 18 full-words for all subprograms. 
(The minimum save area requirements are given in Tables 2 through 6 
for the mathematical subprograms and in Table 16 for the service sub- 
programs. ) 

Aligns the argument fist at a full-word boundary. 

Indicates the first byte of the only argument. 

Contains the address of the argument. 

Contains the address of the first argument. 
Contains the address of the second argument. 



DC 
DC 



X'80' 
AL3(arg„) 



Indicates the first byte of the last argument. 
Contains the address of the last argument. 



Figure 3. General Assembler Language Calling Sequence 



When the load module is executed, the ihcsqrt sub- 
program is called to compute the square root of the 
number in amnt; the result is stored in floating-point 
register 0. The binary calling sequence identifier is 
not used. 

The assembler language user of the Model 44 Pro- 
gramming System will use the calling sequence given 
in Figure 3. He will, however, use only the A-type 
address constant where the choice between that and 
the V-type is given. As the note in the figure states, 
the label in the operand portion of the address constant 
must be made the object of an extrn statement to 
obtain proper linkage. 

Mathematical Subprograms 

The assembler language programmer supplies one or 
more arguments for each mathematical subprogram. 
The arguments may be either integer values or nor- 
malized floating-point real or complex values. 

An integer argument is four bytes in length and 
starts on a full-word boundary. A real argument is 
either four or eight bytes in length. The four-byte 



argumeiit starts on a full-word boundary. The eight- 
byte argument starts on a double-word boundary and 
occupies two adjacent words. The first word contains 
the most significant digits. This word is also the ad- 
dress of the entire argument; the second word con- 
tains the least significant digits. 

A complex argument is either eight or sixteen bytes 
in length and starts on a double-word boundary. The 
first half of the argument contains the real part of the 
complex argument; the second half contains the 
imaginary part. The address of the real part of the 
argument is the address of the entire argument. 

Each mathematical subprogram returns a single 
answer. This answer is either an integer value or a 
normalized floating-point real or complex value. An 
integer answer is stored in general register 0, a real 
answer is stored in floating-point register 0, and a 
complex answer is stored in floating-point registers 
and 2. 

Tables 2 through 6 contain additional information for 
using the mathematical subprograms in an assembler 
language program. These tables give the floating-point 
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Table 17. Assembler Information for the Service Subprograms 



Subprogram 
Name 


Entry 
Name(s) 


Save Area 
(Full Words) 


xxxFDUMP 


DUMP 
PDUMP 


18 
18 


xxxFDVCH 


DVCHK 


10 


xxxFEXIT 


EXIT 


5 


xxxFOVER 


OVERFL 


10 


xxxFSLIT 


SLITE 
SLITET 


9 
10 



registers that are used by the subprogram and the 
save area required by the subprogram. 

Service Subprograms 

The service subprograms do not use the floating-point 
registers during execution; however, each service sub- 
program requires a save area. The minimum size of 
the save area depends upon the subprogram to be 
used and is given in Table 17. 
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Appendix F. Sample Storage Printouts 



A sample printout is given below for each dump for- 
mat that can be specified for the xxxfdump subpro- 
gram. The printouts are given in the following order: 
hexadecimal, logical *i, logical *4, integer *2, 

INTEGER *4, REAL *4, REAL *8, COMPLEX *8, COMPLEX *16, 

and literal (see Figure 4). Note that the headings on 
the printouts are not generated by the system, but 
were obtained by using format statements. 



CALL PDUMP WITH HEXADECIMAL FORMAT SPECIFIED 

00A3E0 "tSSFSElO 00000000 ifSSFSElO 10000000 '(2100000 

006DC8 '♦2BOO0O0 00000000 00000000 00000000 00000000 00000000 OOOOOnOO 00000000 00000000 00000000 00000000 00000000 

a06DF8 COOOOOOO 00000000 1*1200000 1*1566666 OOOOOOOC itUOOOOO 



CALL PDUMP WITH L0GICAL*1 FORMAT SPECIFIED 
006E1E T F 



CALL PDUMP WITH LOGICAL*!* FORMAT SPECIFIED 
006E10 F T 



CALL PDUMP WITH INiTEGER*2 FOR.MAT SPECIFIED 
006E18 10 

006E1A -100 
006E1C 10 



CALL PDUMP WITH INTEGER*!* FORMAT SPECIFIED 

006E20 1 2 

006Ei*8 11 12 



CALL PDUMP WITH REAL*!* FORMAT SPECIFIED 
006E00 0.20000000E 01 0.55999996E 01 



CALL PDUMP WITH REAL*8 FORMAT SPECIFIED 
006DC8 0.1759999999999999D 03 



CALL PDUMP WITH C0MPLEX*8 FORMAT SPECIFIED 

O06DD0 (3.0000000,1*. 0000000) C* . 0000000, 8 . 0000000 ) 



CALL PDUMP WITH C0MPLEX*16 FORMAT SPECIFIED 

006DEO CO. 99999 99999999990, 0.9999999999999990) C-0.99999"n°99999990, -0.9999999999999990) 



CALL PDUMP WITH LITERAL FORMAT SnECIFIED 
006E5C THIS ARRAY CONTAINS ALPHAMERIC DATA 



Figure 4. Sample Storage Printouts 
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